AN  INSTIGATION  OF  THE  EARLY  WAKE 
UliiMIND  SUBMERGED  MAT  PLATES  AS  INFLUENCED 

IJY  TRAILING  l-DGE  CONFIGURATION 


EDWIN  J.   PETERSEN,  JR. 
GERALD  D.  SYLVESTER 


LIBRARY 

O.S.  NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CALIFORNIA 


i 
AN  INVESTIGATION  OF  THE  EARLY  WAKE 

BEHIND  SUBMERGED  FLAT  PLATES 

AS  INFLUENCED  BY 

TRAILING  EDGE  CONFIGURATION 

by 

LIEUTENANT  EDWIN  J.  PETERSEN,  JR.,  U.  S.  NAVY 

B.S.,  U.  S.  Naval  Academy 

(19*3) 

and 

LIEUTENANT  GERALD  D.  SYLVESTER,  U.  S.  NAVY 

B.S.,  U.  S.  Naval  Academy 

(1952) 

SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS 

FOR  THE  MASTER  OF  SCIENCE  DEGREE  IN  NAVAL  ARCHITECTURE 

AND  MARINE  ENGINEERING  AND  THE  PROFESSIONAL  DEGREE, 

NAVAL  ENGINEER 
at  the 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
May,  1961 


Llbrnry 

U.  S.  Naval  Postgraduate  School 

Monterey,  California 

11 

AN  INVESTIGATION  OP  THE  EARLY  WAKE  BEHIND  SUBMERGED  PLAT 
PLATES  AS  INFLUENCED  BY  TRAILING  EDGE  CONFIGURATION 

by 

Lieutenant  Edwin  J.  Petersen,  Jr.,  U.  S.  Navy 

and 
Lieutenant  Gerald  D.  Sylvester,  U.  S.  Navy 

Submitted  to  the  Department  of  Naval  Architecture  and  Marine 
Engineering  on  20  May  1961  In  partial  fulfillment  of  the  re- 
quirements for  the  Master  of  Science  Degree  in  Naval  Archi- 
tecture and  Marine  Engineering  and  the  Professional  Degree, 
Naval  Engineer. 

ABSTRACT 

This  report  presents  the  results  of  an  experimental 
Investigation  of  the  behavior  of  the  turbulent  early  wake 
behind  three  flat  plates  of  varying  trailing  edge  config- 
uration submerged  in  water. 

The  circumstances  surrounding  "singing"  propellers  and 
the  influence  of  vortex  formation  upon  this  and  other 
hydroelastlc  vibration  phenomena  are  discussed.  The  sig- 
nificance of  the  vortex  as  a  forcing  function  is  emphasized. 

The  influence  of  plate  vibration  upon  the  early  wake  is 
essentially  eliminated  in  the  investigation  by  rigidly 
clamping  the  test  plates. 

A  prediction  of  plate  profile  drag  coefficient  and  of 
dimenslonless  shedding  frequency,  the  Strouhal  number,  is 
made  based  upon  von  Karman's  vortex  street  analysis. 

The  longitudinal  and  transverse  distribution  of  local 
mean  velocity,  turbulence  intensity,  and  static  pressure 
were  measured  in  the  early  wake  of  three  selected  test 
plates  of  varying  trailing  edge  configuration.   Turbulence 
energy  density  spectra  were  obtained  at  selected  positions 
in  the  early  wake.  Profile  drag  coefficients  were  computed 
and  vortex  shedding  frequencies  measured,  and  are  compared 
with  predicted  values. 
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I.  INTRODUCTION 

A.   General  Statement  of  the  Problem. 

One  of  the  most  vital  problems  In  ship  design  today,  and 
particularly  in  warship  design,  is  the  reduction  of  noise  and 
vibration.   Of  the  many  and  varied  sources  of  noise  and  vibra- 
tion in  ship  structures,  one  of  the  more  significant  is  the 
hydroelastic  vibration  of  propeller  blades  and  control  sur- 
faces.  The  phenomenon  of  "singing"  propellers  has  interested 
workers  in  several  fields  since  the  1930 's,  when  higher  pro- 
peller speeds  and  the  use  of  airfoil  blade  sections  instead 
of  ogival  sections  became  prevalent1-   .   Fundamentally,  blade 
vibration  of  some  sort  has  been  held  responsible  in  the 
literature  for  this  phenomenon,  but  writers  have  differed  in 
their  analyses  of  the  vibration  and  in  their  several  explana- 

T  1  —Pi   1 

tions  for  the  cause  of  this  vibration1  "  J . 

While  propeller  in-flow  variations,  blade  back  pressure 
variations,  and  wheel  unbalance  may  indeed  contribute  to  the 
phenomenon,  most  writers,  however,  agree  that  the  probable 
cause  of  singing  appears  to  be  the  formation  and  shedding  of 
vortices  in  resonance  with  the  natural  frequency  of  the  blade 
and  its  entrained  water.   If  the  vibrations  set  up  by  the  pro- 
cess of  forming  of  eddies  at  the  trailing  edge  of  a  propeller 
blade  are  in  the  audible  range  of  frequency,  a  singing  note 
may  be  emitted L J J .   It  has  been  found  that  for  a  given  pro- 
peller system  singing  occurs  in  a  relatively  narrow  frequency 
band  and  over  a  relatively  narrow  range  of  propeller  rpm. 

While  audible  singing  is  generally  not  encountered  in 


larger  foil  or  plate  structures  such  as  rudders,  hydrofoils, 
submarine  diving  planes  and  other  control  surfaces,  vibration 
and  flutter  are  significant  problems'-   '    .   It  appears 
reasonable  to  consider  that  vortex  induced  vibration  may  be 
at  least  in  part  responsible  for  this. 

It  has  been  found  in  practice  that  modifications  to  the 
geometry  of  propeller  blade  sections  will  alter  the  character- 

ri2i 

istics  of  singing  or  eliminate  it  entirely.   Den  HartogL   J 
and  others  discuss  the  procedure  of  sharpening  the  leading 
edge  of  propeller  blades  and  turbine  blade  runners  to  eliminate 
undesirable  singing.   Propeller  sections  with  sharp  leading 
edges  have,  however,  been  known  to  sing  violently,  and  several 
authors'-  ,^-'yJ  discuss  the  practice  of  blunting  the  trailing 
edges  of  blades  to  alleviate  this  condition.   The  periodic 
transverse  forces  due  to  alternate  shedding  of  vortices  incite 
a  resonant  mode  of  vibration  in  some  part  of  the  propeller 
blade,  and  according  to  Saunders  L-7J  the  remedy  is  to  change 
the  shape  of  the  trailing  edge  in  some  manner  so  that  the  vor- 
tex formation  is  eliminated  or  that  the  vortices  and  circu- 
lation are  stabilized.   One  simple  method  is  to  cut  off  the 
trailing  edge  square,  so. that  the  separation  points  are  "well 
defined  and  fixed  in  position; "  another  is  to  plane  the 
trailing  edge  flat  at  some  angle  to  the  base  chord.  Various 
other  trailing  edge  geometry  modifications  are  discussed  in 
the  literature. 

Since  practical  experience  shows  that  modification  of 
the  trailing  edge  will  alter  the  response  of  a  blade  section, 
and  since  it  is  not  reasonable  that  such  modifications  would 


materially  alter  the  mass  or  the  Internal  damping  of  the 

T71 
system1 ' J ,  it  follows  that  it  is  the  forcing  function,  i.e., 

the  shedding  of  vortices,  which  has  in  some  way  been  modified. 

A  large  body  of  literature  exists  on  the  subject  of  wake 
flows  behind  various  types  of  bodies.   The  majority  of  it, 
however,  deals  with  the  wakes  behind  bodies  of  revolution  and 
bluff  bodies,  and  almost  all  of  it  treats  the  so-called  fully 
developed  wake  well  downstream  of  the  body.   In  order  to  con- 
duct a  study  of  trailing  edge  vortex  phenomena  at  their 
origin,  it  becomes  necessary  to  investigate  the  "early  wake" 
or  the  region  within  about  ten  characteristic  dimensions 
downstream  of  a  foil  or  flat  plate.   Such  a  study  should  pro- 
perly consider  the  distribution  of  the  energy  present  in  the 
early  wake  and  the  variations  of  velocity,  pressure,  and  turb- 
ulence intensity  in  this  region.   Furthermore,  the  plate  bound- 
ary layer  itself  and  the  free  shear  layer  where  vorticity 
becomes  pronounced  may  be  expected  to  influence  these  phenomena. 
These  considerations  have  prompted  an  investigation  by  the 
authors  of  the  early  wake  behind  three  flat  plates  of  varying 
trailing  edge  configuration  at  the  variable  flow  water  tunnel 
of  the  Massachusetts  Institute  of  Technology  Hydrodynamics 
Laboratory. 
B.   Scope  of  the  Investigation, 

As  part  of  a  general  investigation  of  wake  mechanics 
conducted  at  the  Massachusetts  Institute  of  Technology  Hydro- 
dynamics Laboratory  and  sponsored  by  the  David  Taylor  Model 
Basin,  HuvalL  ~JJ   recently  conducted  an  investigation  of  the 
early  wake  behind  a  fixed  flat  plate  having  a  semi-cylindrical 


leading  edge,  a  chord-to-thickness  ratio  of  8.0  and  a  60° 
included  angle  trailing  edge.   This  particular  plate  was  one 
of  a  series  studied  by  Ippen,  Toebes,  and  Eagleson1  ^"J  for 
the  purpose  of  investigating  hydroelastjic  behavior  of  flat 
plates  with  various  trailing  edge  configurations.   The  authors 
chose  three  other  plates  from  this  systematic  series  for  their 
Investigation;  their  dimensions  and  geometries  are  listed 
below  in  Table  I; 

TABLE  I 
DIMENSIONS  AND  GEOMETRIES  OF  TEST  PLATES 

PLATE  LEADING  EDGE  THICKNESS  CHORD  TRAILING  EDGE 

A  semi-cylindrical     1/J+"  2"      £ZZ3  p  =  90° 

B  semi-cylindrical     l/k"  2"       $         ) 

C  semi-cylindrical     I/I4."  2"      { 


These  three  plates  correspond  respectively  to  Plates  12,  11, 
and  15  of  Reference  [1I4.]. 

The  experimental  program  for  this  investigation  provided 
for  the  measurement  of  velocity  and  pressure  distributions 
and  distribution  of  turbulence  intensity  across  the  wake  at 
several  longitudinal  positions  in  the  early  wake  of  the  test 
plates,  and  for  measurement  of  the  turbulence  energy  spectra 
at  several  longitudinal  positions  on  the  wake  centerline  and 
in  the  vicinity  of  the  paths  of  shed  vortices. 

Plates  A  and  B  were  originally  selected  from  consider- 
ation of  vibrational  amplitude  responses  as  a  function  of 
restraint  (Figure  29  of  1JL4.) j  these  particular  plates  repre- 
sented respectively  the  lowest  and  one  of  the  highest 


5 
response  curves  presented.   During  the  course  of  the  Investi- 
gation, it  became  desirable  to  include  an  intermediate  con- 
figuration, and  Plate  C  was  incorporated  into  the  investigation. 
As  in  Huval's  investigation,  the  authors  decided  to  clamp 
the  plate  rigidly  in  the  water  tunnel  test  section  in  order 
to  be  able  to  study  the  wake  energy  with  minimum  contribution 
from  plate  vibrations 0   In  order  to  facilitate  future  compari- 
son with  Huval's  results,  it  was  decided  to  conduct  measure- 
ments at  a  free  stream  velocity,  U0,  in  the  neighborhood  of 
10  feet  per  second . 


II.  REVIEW  OF  PREVIOUS  WORK 

A.   Vortex  Excited  Blade  Vibration 

Of  the  large  number  of  explanations  postulated  for 
singing  propeller  blades  over  the  past  30  years,  the  sugges- 
tion that  vortex  formation  and  shedding  might  be  the 
responsible  mechanism  appears  to  be  one  of  the  most  plausible. 
Experimental  work  conducted  in  the  United  States  and  in  the 
Soviet  Union  in  recent  years'-  *  -^  support  strongly  the  con- 
tention that  singing  occurs  when  the  vortex  shedding  fre- 
quency is  in  resonance  with  one  of  the  natural  modes  of 
vibration  of  a  given  structure. 

Considering  that  at  resonance  singing  frequency  coin- 

rfli 
cides  with  vortex  shedding  frequency,  Gongwer1  J  presented  an 

empirical  modification  of  the  Strouhal  formulation,  suggesting 

that  the  characteristic  trailing  edge  dimension  should  include 

a  portion  of  the  boundary  layer  thickness .   Introducing  this 

factor,  he  obtained  experimentally  from  his  work  with  vanes 

Sb.=  {&'*  0.185  (1) 

ul 
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and  where  Uj  is  the  component  of  the  free  stream  velocity 
normal  to  the  trailing  edge  of  the  body.   Krivtsov  also 
illustrates'-  'l   the  essentially  constant  behavior  of  Strouhal 
number  for  a  given  trailing  edge  thickness „ 

Several  writers  have  been  concerned  with  the  fact  that 
the  singing  behavior  of  identical  propellers  may  be  quite 


different.   However,  when  one  considers  the  processes  in- 
volved in  propeller  manufacture  and  balancing,  it  is  not  at 
all  surprising  that  two  propellers  of  the  same  class  might 
have  slightly  different  natural  frequencies,  and  that  the 

frequencies  of  the  individual  blades  might  vary  somewhat. 

f7l 
Lankester1,  ■ J  points  out  that  slight  differences  in  natural 

frequencies  between  blades  are  not  at  all  uncommon.  This  may 
increase  the  response  or  decrease  the  equivalent  damping  of 
one  propeller  with  respect  to  another  having  blades  of  iden- 
tical natural  frequency.   If  propeller  singing  is  due  to 
self-excited  vibration,  a  small  difference  in  response  could 
make  all  the  difference  between  singing  and  non-singing. 

Propeller  singing  has  been  shown  to  be  dependent  upon 

rn 

water  temperature1,  J  and  ambient  pressure,  indicating  further 
the  flow  dependence  of  the  phenomenon. 
B.  Vortex  Formation  and  Wake  Periodicity 

Since  Strouhal's  work  in  the  latter  part  of  the  19th 
century,  a  large  number  of  investigators  have  studied  the 
vortex  shedding  of  cylindrical  and  spherical  bodies,  and 
more  recently  of  bluff  bodies.  While  this  work  is  of  con- 
siderable interest  in  theory  and  in  practice,  and  valuable  as 
background,  it  is  not  immediately  applicable  to  the  case  of 
thin  flat  plates  and  other  blade  sections,  and  will  be  treated 
only  cursorily. 

The  Reynolds  number  dependence  of  the  Strouhal  number, 

1*161 

suggested  by  Raleigh1   J  has  been  clearly  demonstrated 

Tl7l 
recently  for  circular  cylinders  by  Kovasznay1  /J  and  by 

r  i  ft  i 

RoshkoL   J.   Roshko  presents  empirical  expressions  for  S^ 
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for  two  Reynolds  number  regimes  of  flows 

Sd  =  0.212  (1  -  21.2/Red)     J+0<Red<l50     (3) 
and 

Sd  =  0.212  (1  -  12.7/Red)    300<Red<£000    (k) 

For  wakes  behind  various  bluff  bodies,  Ippen  e_t  al..  re- 
view the  work  of  several  investigators'-  ^J  and  indicate  that 
the  general  range  of  Strouhal  numbers  based  upon  the  largest 
characteristic  dimension  normal  to  the  direction  of  the  flow 
is  from  about  0.12  to  0.21,  i.e.,  below  the  value  of  0o2 
applicable  to  circular  cylinders. 

Page  and  Johansen  in  experiments  on  the  early  wake  be- 
hind inclined  flat  plates  showed  that  Strouhal  number  is 

dependent  upon  angle  of  incidence,  a,  increasing  as  a  becomes 

flQl 
small L  -/J.   Expressing  S  in  terms  of  the  separation  between 

free  shear  layers,  S  was  found  to  be  approximately  constant 

and  =s  0.16.   Both  in  Gongwer's  experiments  and  in  other  ex- 

periments  on  oval  struts'- ' J  it  was  found  that  the  effect  of 

small  angles  of  incidence  was  to  alter  the  amplitude  of 

vibration  without  materially  affecting  the  frequency. 

In  none  of  the  theoretical  treatments  of  the  problem 
is  an  explanation  of  the  actual  mechanism  of  vortex  formation 
presented.  The  work  of  BetchovL   J  suggests  that  a  possible 
approach  may  lie  in  the  treatment  of  boundary  layer  oscil- 
lations. 

Von  Karman's  classical  treatment  of  periodic  wakes  has 
crystallized  and  stimulated  significant  work  in  this  field. 
His  stability  criterion  for  a  vortex  array  in  a  wake  has 
theoretical  and  practical  significance  and  will  subsequently 


be  utilized  to  predict  shedding  frequencies  and  to  analyze 
free  shear  layer  separation.  This  concept  and  Birkhoff's 
modification  will  be  discussed  further  under  "Theoretical 
Considerations,"  Von  Karman's  treatment  makes  no  attempt  to 
describe  the  mechanism  of  vortex  formation  and  is  in  fact 
independent  of  the  mechanism  and  of  the  geometry  of  the 
structure  behind  which  it  applies. 

C.  Wake  Structure 

The  similarity  principles  of  von  Karman  applied  to  the 
Reynolds  equations  of  motion  have  proven  useful  in  treating 
axially  symmetric  or  two-dimensional  turbulent  wakes.  Un- 
fortunately, similarity  hypotheses  apply  to  a  state  of 
statistical  equilibrium  of  flow  processes  relatively  far 
downstream  in  a  wake,  and  cannot  be  considered  applicable  to 
the  early  wake, 

r-iOi  IT  9"l  1 

Roshko's  experiments1,   J  and  those  of  Cooper  and  LutzkyL   J 
present  results  of  velocity  and  turbulence  distributions  in 
the  relatively  early  wake  behind  cylinders  and  bluff  bodies 
respectively, 

r  99 1 

BirkhoffL   J  and  others  point  out  that  contrary  to  the 
predictions  of  the  Kirchoff  wake  model,  actual  large  wake 
underpressures  exist. 

D,  Trailing  Edge  Geometry 

Ippen  e_t  al.  point  out  the  effect  of  trailing  edge 
geometry  on  drag  coefficients  for  various  shapes,   Schlieh- 
tingL  -IJ  shows  that  the  drag  coefficient  can  be  related  to 
the  vortex  street  through  momentum  analysis,  but  notes  that 
this  procedure  is  applicable  only  relatively  far  downstream 
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where  the  static  pressure  at  the  measuring  station  is  equal 
to  the  free  stream  pressure.,  The  correction  methods  of  Betz 
and  of  Jones'- ™""  may  be  employed  in  the  early  wake  region, 
and  have  yielded  valid  results  behind  airfoil  sections  as 
close  to  the  trailing  edge  as  0.05  chord  lengths. 
E.   Turbulence 

The  complex  nature  of  wake  turbulence  and  turbulent 
(or  even  periodic)  wake  energy  precludes  any  detailed  dis- 
cussion here.   Reference  [28]  in  particular  and  portions  of 
references  [29]  through  [31]  are  of  interest  for  this  inves- 

rr^i 

tigation  as  are  the  results  of  Huval's  experiments1,  JJ. 

Ambient  turbulence  background  measurements  were  made 

for  the  water  tunnel  test  section  by  Perkins fcJ  J  and  by 

Fl^l 
Huval"-  •JJ   indicating  relative  mean  turbulence  levels  of 

about  0.03  and  0,065  respectively.  Assuming  isotropic  tur- 
bulence with  a  general  power  decay  law,  and  ergodic  process, 
Huval's  measurements  indicated  that  the  free  stream  turbulence 
was  largely  in  its  initial  stage .  The  character  of  the  free 
stream  turbulence  spectrum  obtained  by  Huval  indicates 
ani sot ropy  in  the  low  frequency  end;  the  high  frequency  range, 
however,  indicates  isotropic  statistical  equilibrium.  Huval 
notes  that  the  high  frequency  spectrum  appears  to  follow  von 
Karman's  suggested  energy  variation  with  yT^'j   "■  ™  .  This 
is  in  agreement  with  the  Kolmogoroff  spectrum  law  cited  by 
Hinze[28]. 

Hinze  points  out  that  in  all  cases  where  the  mean  vel- 
ocity shows  a  gradient  turbulence  will  be  anisotropic »   It 
should  then  be  expected  that  in  regions  of  a  wake  where  such 
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gradients  exist  anisotropic  turbulence  will  obtain.   This  con- 
dition is  expected  to  exist  in  the  early  wake  behind  an  object. 

r  -l  d  I 

Roshko1   J  notes  that  the  velocity  fluctuation  in  the 
wake  of  a  vortex- shedding  cylinder  displays  a  predominant 
frequency  which  is  the  shedding  frequency]  he  points  out 
further  that  this  predominant  frequency  will  have  random 
irregularities  imposed  upon  it. 

Distributions  of  turbulence  intensity,  mean  velocity, 
and  pressure  in  a  turbulent  wake  are  discussed  by  Hinze. 
Theoretical  treatments  of  these  distributions  are  in  general 
predicated  upon  similarity  assumptions  which  are  well  known. 
Such  assumptions,  however,  are  not  expected  to  be  Justified 
in  the  early  wake  behind  a  structure.   Cooper  and  Lutzky11   J 
found  that  certain  predictions  made  on  the  basis  of  similarity 
hypotheses  were  not  upheld  by  experimental  evidence  in  a 
region  relatively  close  behind  bluff  bodies.   They  note  that 
this  discrepancy  was  not  resolved. 
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III.  THEORETICAL  CONSIDERATIONS 

A.   Analysis  of  the  Vortex  Street, 

From  consideration  of  complex  potential  theory,  It  can 
be  shown1  ^J  that  the  complex  potential  function  due  to  an 
infinite  row  of  equally  spaced  vortices  is  given  by 

w  =  %  in  sin  22.  (5) 

from  which  it  can  further  be  shown  that  for  two  such  rows 
parallel  to  each  other  and  separated  by  a  distance  b, 

w  =  ^  [in  sin  £  (z  -  |)  -  In  sin  £  (z  -  f  +  ^)]  (6) 

where i 

a  =  longitudinal  vortex  spacing 
b  =  transverse  vortex  spacing 
k  =  vortex  strength 

dw 
Since  u  -  iv  =  -4=-  ,  differentiation  and  separation  of 

components  yields  an  expression  for  vortex  velocity  in  the 

longitudinal  directions 

iv  =  £  tanh  - 


uv  =  £  tanh  &  (?) 


Prom  stability  considerations,  von  Karman  showecr  ^'   ■* 
that  the  street  can  be  stable  only  when  vortices  in  one  row 
are  opposite  the  interval  center  between  vortices  in  the  other 
row,  and  then  only  when  the  ratio  b/a  =  0.281. 

Rosenhead  (cited  in  25)  shows  that  rigid  boundaries 
equidistant  from  the  mean  centerline  of  the  vortex  street 
may  affect  the  stability  ratio.   If  h  -  the  distance  between 
the  rigid  boundaries,  then  as  the  ratio  a/h  varies  from  0  to 
0.8l£,  the  stability  ratio  has  a  unique  but  continuously 
changing  value,  going  from  0.281  to  0.256.  For  values  of  a/h  > 
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0.81^  stability  exists  for  a  certain  range  of  values  of  b/a, 

and  for  values  of  a/h  >  l.i+19.  the  configuration  is  stable 

for  all  values  of  b/a.   Birkhoff  (cited  in  llj.)  predicts  a 

stability  ratio  of  b/a  =  0.35  from  inertia  considerations. 

The  test  section  of  the  water  tunnel  used  in  this  present 
investigation  is  rectangular,  measuring  9  inches  high  by  7-1/2 
inches  wide.  The  test  plate  is  mounted  vertically  in  the 
section.   The  dimension  h  between  rigid  boundaries  in  this 
case  is  then  7-1/2".   Assuming  that  at  the  trailing  edge 
b  =  plate  thickness,  t,  and  proceeding  from  the  von  Karman 
criterion  a/b  =  0.281  as  a  starting  point,  the  ratio  a/h  = 
0.1185.  Following  a  reiterative  procedure  and  assuming 
linear  decrease  of  stability  ratio  with  increase  in  a/h,  the 
Rosenhead  treatment  indicates  that  the  stability  ratio  for 
vortex  streets  in  the  current  investigation  would  be  b/a  = 
0.271.   In  actuality  the  transverse  spacing  is  expected  to 
be  greater  than  the  plate  thickness  due  to  wake  spreading. 
The  assumption  of  linear  decrease  in  b/a  is  considered 
adequate  in  the  region  of  a/h  concerned.   It  appears,  then, 
that  a  stability  ratio  for  b/a  may  lie  in  the  region  from 
about  0.27  to  about  0.28  for  this  investigation. 

After  the  manner  of  Ippen^1^-"  and  Birkhoff  °-27-"  it  can 
be  shown  that  considerations  of  vorticity  diffusion  from  the 
boundary  layer  can  lead  to  an  expression  for  the  strength,  k, 
of  the  vortices  in  the  wake.   If  K^  be  the  rate  of  shedding 
vorticity  from  each  side  of  an  obstacle,  then  as  an  approx- 
imations 

Ki  =   J'cuxdy  =  j  |H.  Uldy  =  ^i*  (8) 
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where  \i^  is  the  velocity  outside  the  boundary  layer  near  the 
separation  point.  If  £  =  K/K]_,  where  K  is  von  Karman's  rate 
of  vorticity  discharge  from  each  side  of  an  obstacle,  then  /3 
is  a  measure  of  the  anihllation  of  vorticity  at  the  trailing 
edge.  Both  Prandtl  and  Roshko( cited  in  li^.  and  27)  have  found 
experimentally  that  £  =  0.5  for  periodic  wakes.  Applying  this 
value  and  assuming  that  ui  =  U,  we  obtains 

K  =  0.25  U2  (9) 

Introducing  this  into  von  Karman's  expression  for  vortex 
strength,  we  obtains 

Introducing  this  into  Eq.  (7)  we  obtains 


u. 


gZ  (1   -  jjZ)   =  0.125  tanh  trb/a  (ll) 

from  which: 

^  =  0.5(1   -Vi  -  0.5  tanh  Trb/a   )  (12) 

Taking  as  the  stability  ratio  b/a  =  0.271  and  operating 
on  Eqs.  (10)  and  (ll)  we  obtains 

uv/U  =  0.0955  (13) 

k  =  0,277  aU 
From  von  Karman's  relationship  between  profile  drag,  Dp, 
and  wake  vorticity,  where  Dp  =  kp|-  (U-2uv)  +  £]■—■,  and  from 
Eqs.  (13)  we  obtains 

Dp  =  0o0728  paU2  (li|.) 

from  which: 

cDp  =  fpTO^  °»1456  aA  (15) 
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and  reintroducing  the  stability  ratio  b/a  =  0.271s 

CDp  "  °»537  b/t  (16) 

It  should  be  pointed  out  that  Ippen  e_t  al.  assumed  a 
value  of  K/u2  «  0.I4.  which  is  in  consonance  with  empirical 
results  cited  by  Birkhoff  *■.   If,  however,  the  experimental 
result  of  jS  =  0.5  is  accepted  and  Eq.  (8)  is  considered  valid, 

this  would  lead  to  a  value  of  u^/D*  =  1.265? • 

r  271 

Birkhoff1  ' J  points  out  that  it  is  a  corollary  con- 
sequence of  the  theorem  of  conservation  of  momentum  that  in 
viscous  wakes  vortex  transverse  spacing  tends  to  increase; 
this  is  experimentally  observed.   For  slightly  viscous  fluids 
it  is  further  noted  that  transverse  vortex  spacing  is  essen- 
tially equal  to  the  wake  thickness  a  little  behind  the  ob- 
stacle.  Birkhoff  notes  that  for  circular  cylinders  the 
critical  diameter  Reynolds  number  above  which  the  wake  first 
becomes  periodic  is  lowered  by  external  turbulence.  For  flat 
plates  parallel  to  the  flow  the  critical  wake-thickness  Rey- 
nolds number  is  about  3=W4-\/^el  where  Rei  =  "v~  »   this  is  the 
same  order  of  magnitude  as  in  the  case  of  a  circular  cylinder. 

Following  Ippen' s  procedure,  the  plate  thickness  Strouhal 

number  for  vortex  shedding  can  be  estimated  from  Eqs.  (13) 

and  (15): 

«         nt  _   M        uvw,/al  _  0.i3t2 


Introducing  Eq.    (16)   we  obtain: 

t 

b 


St  =  0.2^5  I  (17) 
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from  which: 

StU   0.2J4-5  U  nftN 

nv  -  -g ^—  (18) 

If  wake  thickness  and  transverse  vortex  separation  are  taken 
as  being  approximately  equal  to  plate  thickness  in  the  early 
wake,  then  for  t  =  0.25"  w©  obtain  a  predicted  vortex  shed- 
ding frequency  of; 

nv  =  11.76  U  (cps)  (19) 

It  should  be  pointed  out  that  none  of  the  foregoing 
takes  into  consideration  trailing  edge  geometry.   Furthermore, 
these  developments  do  not  consider  any  Reynolds  number  de- 
pendence upon  either  Strouhal  number  formulation  or  periodic 
wake  structure. 
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IV.  DESCRIPTION  OF  EQUIPMENT  AND  INSTRUMENTATION 

A.  The  Variable  Flow  Water  Tunnel. 

A  detailed  description  of  the  variable  flow  water  tunnel 
is  given  in  Reference  [33] •   Photographs  of  the  tunnel  test 
section  are  presented  as  Figs.  la.  and  lb.;  Fig.  Ic.  is  a 
schematic  diagram  of  the  tunnel  and  its  associated  flow  cir- 
cuit.  The  hydraulic  characteristics  of  the  water  tunnel  are 
presented  in  Appendix  A  of  Reference  [li|.]a   Tunnel  operation 
is  possible  with  any  one  or  a  combination  of  three  centri- 
fugal pumps. 

Considerations  of  minimum  mean  level  of  background 
turbulence  and  consistency  of  pump  operation  led  the  authors 
to  conduct  all  measurements  with  Pump  III  providing  the  flow. 
Referring  to  Fig.  Ic,  the  flow  was  routed  from  Pump  III  to 
the  right  through  the  18"  distribution  circuit,  back  to  the 
left  past  Pump  I  (which  was  isolated  from  the  circuit),  and 
up  through  the  water  tunnel  inlet  riser.   It  had  initially 
been  planned  to  operate  with  both  Pumps  II  and  III  in  oper- 
ation simultaneously,  with  the  stop  valve  between  them  shut, 
but  flow  reversal  through  Pump  II  and  problems  in  velocity 
control  led  to  the  decision  to  isolate  Pump  II  from  the 
circuit. 

From  the  expanding  tunnel  inlet  riser  the  flow  passes 
through  four  No.  16  mesh  screens  and  a  33"  long  contraction 
into  the  test  section.   The  plexiglass  test  section  is  36" 
long  and  rectangular  in  section,  measuring  9"  high  and  7-1/2" 
wide.   From  the  test  section,  the  flow  returns  to  the 
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FIGURE  la.  Variable  Plow  Water  Tunnel 
Test  Section,  Showing  Location  of  Amp- 
lifier and  Preamplifier. 


FIGURE  lb.  Variable  Flow  Water  Tunnel 
Test  Section,  Left  Side. 
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storage  tank. 

B.  Test  Plates. 

The  test  plates  studied  In  this  investigation  are 
described  briefly  in  Table  I  and  in  more  detail  in  Reference 
[1I4.].   They  were  manufactured  from  1/V  aluminum  stock  and 
carefully  machined.  The  plates  spanned  the  9"  vertical 
dimension  of  the  test  section,  and  presented  a  2"  chord  in 
the  direction  of  the  flow.   The  test  plates  were  anodized  to 
prevent  corrosion  and  pitting!  this  process  proved  to  be 
quite  successful. 

C.  Test  Plate  Mounting. 

The  test  plates  were  mounted  using  the  spindles  illus- 
trated in  Figure  6  of  Reference  [1I4.].   The  spindles  were 
attached  to  the  plates  with  set  screws  such  that  the  center- 
line  of  the  spindles  was  essentially  in  line  with  the 
center  of  radius  of  the  semi-cylindrical  plate  leading  edge. 
The  spindles  were  clamped  to  prevent  rotation. 

Each  of  the  spindles  passed  through  a  roughly  rectan- 
gular 1"  by  2-1/2"  hole  in  the  upper  and  lower  test  section 
walls.  These  holes  opened  further  into  larger  cavities.  The 
upper  cavity  was  originally  sealed  with  a  solid  cover  platej 
because  of  the  tendency  of  air  bubbles  to  accumulate  in  this 
cavity,  a  vent  was  subsequently  installed  in  the  cover  plate. 
The  lower  spindle  passed  through  a  gasketed  circular  hole  in 
the  bottom  cover  plate  and  was  again  braced  below  the  test 
section ,   This  arrangement  is  due  to  the  earlier  investigation 
reported  in  Reference  [ll+]. 


21 

D.  Traversing  Mechanism. 

In  order  to  obtain  data  on  transverse  and  longitudinal 
wake  structures  with  the  flow  measurement  devices  to  be 
described  subsequently,  it  was  necessary  to  provide  for 

traversing  in  these  directions.   The  traversing  mechanism  used 

fi-31 

by  Huval1  JJ  was  employed  in  this  investigation.   The  in- 
strumentation strut,  located  just  downstream  of  the  test 
plate,  passes  through  the  test  section  side  walls  through 
slotted  blocks  which  permit  a  longitudinal  traverse.   These 
are  shown  installed  in  Figs.  Ila  and  lib.  The  instrumentation 
strut  itself,  illustrated  in  Fig.  lie,  can  be  traversed  in 
the  transverse  direction  through  the  side  plates  and  blocks. 

Longitudinal  positioning  of  the  strut  is  accomplished 
by  means  of  scales  graduated  in  I/6I4.  inch  mounted  on  the 
slotted  blocks.  A  transverse  scale  graduated  in  0.002  inch 
permits  transverse  positioning. 

E.  Flow  Measurement  Devices. 

The  total  and  static  head  tubes  and  the  turbulence 
pressure  transducer  employed  in  the  wake  surveys  are  described 
in  detail  in  References  [13]  and  [32].  These  instruments 
permit  measurement  of  the  total  and  static  pressures  in  the 
wake  behind  the  test  plate,  and  will  be  treated  separately 
below.   In  addition  to  these  devices,  a  contraction  manometer 
provides  for  measurement  of  mean  free  stream  velocity. 

1.  Turbulence  Pressure  Transducer. 

A  cylindrical  disc-shaped  barium  titanate  ceramic 
transducer,  0.125"  in  diameter  and  0.100"  thick,  is  mounted 
in  a  tapered  water-filled  probe  tip  illustrated  in  Fig.  III. 
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FIGURE  Ha.  Left  Side  View  of 
Test  Section,  Showing  Position 
of  Test  Plate  and  Transverse 
Strut  Positioning  Scale. 


FIGURE  lib.  Right  Side  View  of 
Test  Section,  Showing  Total  and 
Static  Head  Leads  and  Transdu- 
cer Cable  to  Preamplifier. 
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FIGURE  lie.  Photographs  of  Traversing  System. 
(Prom  Huval,  "The  Early  Wake  of  a  Fixed  Flat 
Plate.") 
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The  ceramic  is  located  behind  a  thin  silver-plated  brass 
diaphragm.   The  tapered  probe  tip  ends  in  a  1/16"  diameter 
with  a  0.025"  hole  drilled  in  it  to  permit  "point"  pressure 
fluctuation  measurements.   Pressure  variations  are  transmitted 
to  the  ceramic  through  a  water  column  expanding  from  0.025" 
diameter  to  approximately  0.221"  diameter,  and  through  the 
diaphragm. 

Considerable  care  must  be  taken  to  insure  that  no  air 
bubbles  are  present  in  the  water  column.   Experience  has  shown 
that  the  presence  of  even  minute  air  bubbles  will  materially 
affect  the  response  of  the  system  by  altering  the  response 
of  the  water  column.   To  this  end,  deaerated  water  is  used  for 
filling  the  probe  tip.   Reasonable  success  was  experienced 
by  pumping  deaerated  water  through  the  probe  tip  for  several 
minutes  with  a  syringe. 

It  was  found  earlier  that  the  transducer  did  not  respond 
to  static  pressure  changes  L *\   o   Recent  improvements,  however, 
have  improved  the  response  of  the  transducer  and  made  quasi- 
static  calibration  possible.   This  will  be  discussed  further 
in  Appendix  A. 

2.  Transducer  Response. 

The  Instantaneous  response  of  a  total  head  tube  aligned 
parallel  to  the  direction  of  mean  flow,  in  Huval's  notation  ^ ^ 
is  given  by? 

Htg  =  p  +  pj_  +  (u  +  u'  )2 

Y  Y       2g 

=  P  +  P?  +  u2  +  2uu°  +  u° 2  (20) 

Y  Y     2g     2g     2g~" 

assuming  that  the  terms  involving  vf  and  w*  do  not  significantly 
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contribute  to  the  total  instantaneous  head  measured  in  the 
longitudinal  direction. 

Now  let  etK  =  the  instantaneous  output  signal,  of  the 
turbulence  gage  in  millivolts,  and 
Kt_  =  turbulence  gage  sensitivity  in  millivolts 
per  foot  of  water. 
Then  the  turbulence  gage  output  signal  becomes: 

etg  =  KtgHtg  =  Ktg(/p  +  pi  +  E  +  ^ii  +  slif)         (21  ) 

\Y    Y     2g    g     2g  / 
Now  let  Gp  =  preamplifier  gain,  and 

e  =  voltage  output  of  the  preamplifier  in  milli- 
volts, 
Then: 

ep  =  GpEtg  f?   +  p_^  +  u^  +  uu^_  +   u'2]  (22) 

\Y   Y    2g    g    2g  / 

furthermore : 

ep  =  ep  +  ep' 

•p  =  avKts(l  +  E!  +  EI^ 

\Y    2g    2g  J 

Referring  to  the  electronics  equipment  block  diagram 

(Pig.  V)  it  is  noted  that  a  DO  filter  (i.e.,  a  1  MFD  capacitor) 

is  installed  between  the  preamplifier  and  the  amplifier,  since 

only  the  fluctuating  components  of  the  turbulence  gage  signal 

are  of  interest  at  the  measuring  instruments.   Therefore,  the 

input  signal  to  the  amplifier  becomes: 

e  '  -  e   -  e 
eP   ~  eP    eP 

=  GpKtg  /pj_  +  uu7_  +  u'2  -  u"7^  \ 

\T"     S       2g     )  (23) 

Letting  em  =  amplified  turbulence  signal  input  to  the  measur- 
ing Instruments,  and 
Ga  =  amplifier  gain,  then 


26 


em  =  GaGpKtg  /V  +  uuj_  +  u'2  -  "'2 


(23) 


\Y      g        2g     / 
Here  em  represents  the  amplified  turbulence  gage  signal  in 
millivolts  at  the  input  to  the  measuring  Instruments  (rms 
voltmeters  and  vibration  analyzer) . 

The  true  rms  voltmeter  output  (measured  by  the  thermo- 
couple voltmeter)  becomes: 


rms  =vem2  =  &asp*tg 


— 

w 

+  uu' 
~g~ 

+  u'2  -  jri\2 

2g          / 

*     (24) 


Expansion  of  the  terms  within  parentheses  in  Eq„ (24)  results 
in: 

A>»  +  uu'  +  u"2  -  u^V  =  u2u'2  +  uu'(u'2  -  u"1"2) 


2g 


g' 


+  (u'2  -  u,2)2  +  2uu'p 


4g2         Yg 
+  (u'2  -  u^p'  +  pj_2 
Yg  Y2 

The  full  expression  for  rms  voltmeter  output  becomes: 


vrms  =  G«SrEt.gf  ugu12  +  uu°3  +  uj^  -  u'2  u«2 


g' 


2g« 


2g: 


+  2uu'p'  -  u'2p°  -  p°2\s 

Y2  / 


(25) 


Yg       Yg 

Now  if  it  may  reasonably  be  assumed  that  u'«u,  say  for 
example  that  V^?   ^.0.1  u,  and  that  p'/Y  is  of  the  order 

of  u1  .  only  the  first  term  within  the  parentheses  makes  a 

2g 
significant  contribution  to  the  rms  signal.   Hence: 

vrms  *  GaGpKtg  Sn^T  (26) 

g 
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3.  Total  Head  and  Static  Pressure  Tubes. 

The  total  head  tube  employed  Is  In  fact  a  standard  Pitot- 
static  tube,  but  only  the  dynamic  tap  Is  utilized.   Because  of 
the  longitudinal  pressure  gradients  present  in  the  wake  of  the 
test  plate  and  their  effect  upon  measurements,  the  static  taps 
of  the  total  head  tube  could  not  be  used,  and  a  separate  static 
head  tube  was  employed.   The  pressure  taps  of  the  static  head 
tube  were  in  line  with  the  tips  of  the  total  head  tube  and  the 
transducer  probe.   See  Flg„  III.   Huval L  <J  gives  a  detailed 
description  of  these  instruments.  Both  pressure  tubes  are 
connected  to  mercury  manometers. 

4.  Contraction  Manometer. 

Connected  to  the  upstream  and  downstream  ends  of  the  water 
tunnel  contraction  are  a  mercury  manometer  and  a  manometer 
filled  with  blue-colored  Merlam  No.  D-8325  Indicating  fluid 
(specific  gravity  =  1.75)°   These  manometers  and  their  piezo- 
meter taps  are  described  in  Appendix  A  of  Reference  U^J  .   The 
contraction  manometer  was  used  to  duplicate  flow  conditions  in 
the  water  tunnel  on  separate  days  of  operation.   Calibration  of 
this  manometer  was  necessary  for  the  particular  flow  conditions 
experienced,  and  is  discussed  subsequently. 
F.  Electronics  Equipment. 

Electronics  equipment  utilized  in  the  measurements  consis- 
ted of  an  electrometer  preamplifier,  a  KinTel  Model  1 1 1  BF 
amplifier,  a  General  Radio  Type  762-B  vibration  analyzer^  a 
Ballantine  vacuum  tube  voltmeter  and  a  Sensitive  Research  Instr, 
Corp.  Model  A  thermocouple  voltmeter,  and  a  Dumont  Type  333 
dual-beam  cathode  ray  oscillograph.   A  Hewlett-Packard  Model 
200C  audio  oscillator  was  employed  in  equipment  calibration. 
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PIGURE  III.  Details  of  Turbulence  Probe  and  Pressure 
Tubes.  (Prom  Huval,  "The  Early  Wake  of  a  Plied  Plat 
Plate.") 
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These  components  are  illustrated  in  Figs.  IVa  and  IVb, 

(At  the  end  of  the  investigation  a  General  Radio  Model 
1445_A  vibration  analyzer  became  available;  this  instrument 
extended  the  frequency  range  of  spectrum  measurements  and 
was  used  in  one  location  in  the  wake  of  Plate  C  for  comparison 
purposes, ) 

It  had  originally  been  intended  to  record  the  amplified 
signal  from  the  pressure  transducer  on  magnetic  instrumentation 
tape  for  remote  spectrum  analysis  on  another  wave  analyzer,, 
Recording  difficulties  made  this  impossible  during  the  course 
of  this  investigation,, 

A  block  diagram  of  the  electronics  components  is  pre- 
sented as  Pig,  V. 

1 ,  Electrometer  Preamplifier, 

The  primary  purpose  of  the  electrometer  preamplifier  is 
to  provide  high  input  impedance  to  the  first  stage  of  ampli- 
fication.  Huval  L1^  describes  the  circuitry  of  this  component, 
and  points  out  that  its  addition  to  the  system  should  permit 
response  to  frequencies  below  1  cps„   The  measured  gain  of  this 
component  (see  Appendix  A)  is  0,4-  over  the  majority  of  the 
operating  frequency  range, 

2,  Amplifier, 

The  EinTel  amplifier  used  in  the  system  has  a  very  flat 
response  over  the  operating  range  of  frequencies.   The  design 
gain  accuracy  with  a  fixed  gain  setting  is  *  \%   from  DC  to  2  ECS, 
The  measured  response  characteristics  (see  Appendix  A)  over  the 
operating  frequency  range  point  this  out  further,   A  DC  filter 
is  connected  to  the  input  of  the  amplifier. 
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FIGURE  IVa.  Electronics  Equip- 
ment: Analysis  Components. 


FIGURE  IVb.  Electronics  Equip- 
ment: Preamplifier  and  Ampli- 
fier Installation. 
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FIGURE  V  Electronics  Equipment  Component  Block  Diagram 
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3.  tfave  Analyzer. 

The  General  Radio  Type  762-B  vibration  analyzer  is  a 
constant  percentage  band  width  tunable  voltmeter  with  a 
logarithmic  scale.   Its  operating  range  is  from  2.5  to  750 
cps.   For  the  sharp  selectivity  position  (in  which  all  measure- 
ments were  made)  the  relative  attenuation  is  3  db  at  a  frequency 
differing  by  \%   from  tuned  frequency.   Calibration  was  per- 
formed using  a  signal  from  a  Hewlett-Packard  audio  oscilla- 
tor, and  is  discussed  further  in  Appendix  A.   A  damping 
capacitance  was  inserted  in  parallel  with  the  analyzer  meter 
to  reduce  needle  fluctuations. 

The  Model  1 554- A  analyzer  used  briefly  at  the  end  of  the 
investigation  is  a  similar  instrument,  but  has  an  extended 
high  limit  of  25,000  cps  for  its  operating  frequency  range. 

4.  Voltmeters, 

The  Ballantine  VTVM  is  a  true  rms-indicatlng  meter  only 
for  sinusoidal  waveforms.   The  thermocouple  voltmeter,  on  the 
other  hand,  will  indicate  true  rms  voltages  for  any  signal 
input . 

5.  Oscilloscope. 

The  Dumont  dual-beam  oscillograph  was  used  primarily  for 
visual  monitoring  of  the  amplified  turbulence  signal.   In  this 
respect  it  proved  very  useful  in  determining  the  presence  of 
air  bubbles  in  the  turbulence  probe  tip.   Such  bubbles  attenuate 
severely  the  higher  frequencies  of  pressure  fluctuation.   See 
Pigs.  Via,  VIb,  and  Vic. 
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FIGURE  Via.  Turbulence  Signal  FIGURE  VIb.  Turbulence  Signal 
Presentation  in  Wake  of  Plate  Presentation  on  Wake  Center- 
ed, y/t=0.5»  Illustrating  Pre-  line,  Illustrating  Predomln- 
dominant  Vortex  Frequency.  ant  Double  Shedding  Frequency. 


FIGURE  Vic.  Turbulence  Signal 
Presentation  with  Air  Bubble 
in  Transducer  Probe  Tip,  Same 
Location  as  Figure  VIb  above. 
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V.  CONDUCT  OF  THE  INVESTIGATION 

A.  Ambient  Turbulence. 

A  study  was  made  of  the  ambient  mean  level  of  turbulence 
as  indicated  by  the  amplified  rms  voltage  from  the  pressure 
transducer  for  various  modes  of  pump  operation.   As  noted 
earlier,  Pump  III  and  the  previously  described  flow  circuit 
were  used  in  the  final  investigation;  this  combination 
yielded  one  of  the  lowest  mean  levels  of  turbulence. 

It  had  been  decided  to  operate  the  tunnel  at  a  flow 

velocity  of  approximately  10  fps  to  facilitate  correlation 

M-ai 
with  Huval's  results1  JJ.   It  subsequently  developed  after 

final  calibration  of  the  contraction  manometer  that  measure- 
ments had  been  taken  at  a  value  of  free  stream  velocity  in 
the  test  section,  U0,  of  about  9*5  fps. 

Energy  spectrum  measurements  were  made  for  the  empty 
test  section  at  one  longitudinal  location  on  the  tunnel 
centerline. 

B.  Wake  Survey  Measurements. 

Pressure,  mean  local  velocity,  and  turbulence  intensity 
distributions  across  the  wake  were  desired  at  representative 
locations  in  the  wakes  of  the  three  test  plates,   A  coordinate 
system  was  defined  as  follows s   the  x-axis  was  taken  as  the 
flow  centerline,  which  was  assumed  to  coincide  with  the  plate 
and  tunnel  centerlines.  The  +x-direction  was  taken  in  the 
downstream  direction,  measured  from  an  origin  located  in  the 
transverse  plane  of  the  plate  trailing  edge.   The  origin  of 
the  y-axis  was  taken  at  the  particular  longitudinal  station 
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where  measurements  were  desired  (x/t  =  I.50,  3. 25.  6.50.  and 
9.75).  The  +y-direction  was  taken  in  the  transverse  (hori- 
zontal) direction  toward  the  right  side  of  the  test  section 
(see  Fig.  la.).   It  was  decided  to  measure  total  and  static 
head  and  turbulence  rms  voltage  across  the  wake  in  increments 
of  0.1  plate  thickness.  Full  wake  surveys  were  made  from 
y/t  =  -1.5  to  y/t  =  +1.5  in  this  manner  at  longitudinal 
stations  x/t  =  3. 25  and  9. 75.   Half-wake  surveys  (from  the 
flow  centerline  to  y/t  =  +1.5)  were  made  at  x/t  =  1.50  and 
6.50. 

C  Wake  Energy  Spectrum  Measurements, 

Turbulence  rms  voltage  and  vibration  analyzer  measure- 
ments were  taken  at  y/t  =  0  (flow  centerline)  and  y/t  =  +0.5 
(in  line  with  the  plate  edge)  at  each  of  the  longitudinal 
measuring  stations  mentioned  above,  for  each  of  the  three 
test  plates.   The  empty  test  section  background  spectrum 
mentioned  in  Section  V.A0  above  was  taken  at  y/t  =  0  and 
xA  =  3.25. 
D.  Longitudinal  Variation  Measurement. 

In  addition,  it  was  desired  to  study  the  longitudinal 
variation  of  pressure,  velocity,  and  turbulence  intensity 
along  the  centerline  and  at  y/t  =  +1.5.  To  this  end  the 
measurements  described  in  Section  v.B.  were  also  conducted 
at  these  transverse  locations  at  the  following  longitudinal 
locations,  x/t  =  0.5.  0.75,  1.0.  1.5,  2.0,  2.75.  3.25.  I>.0. 
5.0,  6.5,  7.5,  8.5,  and  9.75. 
_E.  Calibration  of  Equipment. 

Calibration  of  equipment  is  discussed  in  Appendix  A. 
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F.  Flow  Duplication. 

With  the  foot  valve  at  the  discharge  end  of  the  water 
tunnel  set  and  left  in  one  position,  it  was  found  possible 
to  adjust  the  pump  discharge  valve  and  the  distribution 
circuit  stop  valve  to  duplicate  very  closely  both  velocity 
and  pressure  conditions  within  the  test  section.   The  Meriam 
fluid  contraction  manometer  was  used  for  this  purpose.   This 
procedure  was  followed  in  establishing  flow  during  the 
measurements  made  for  Plates  B  and  Cj  minor  changes  were 
made  in  the  foot  valve  setting  during  the  Plate  A  measure- 
ment s . 
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VI.  REDUCTION  OF  DATA 

A.  Mean  Flow  Velocity,  Up, 

The  contraction  manometer  was  used  to  measure  free  stream 
velocity.  The  Meriam  blue  fluid  manometer  lent  itself  to 
more  precise  reading  and  was  thus  used  primarily  rather  than 
the  mercury  manometer.  The  contraction  manometer  was  cali- 
brated against  a  Prandtl-pitot  tube  inserted  through  the  upper 
tunnel  wall  at  the  upstream  end  of  the  test  section.  Velocity 
of  the  free  stream  was  obtained  from  the  expression: 

U0  =  2,06V  AhM  (27) 

where  Ah^  =  Meriam  blue  fluid  manometer  difference 
reading  in  inches. 

B.  Local  Mean  Velocity,  u. 

Local  mean  velocity  was  computed  from  the  total  and 
static  head  tube  manometer  readings  using  the  expressions 


U  ~/2g  Y&  (Aht  -  Ahs)  +  Az 


(28) 


where?    Aht  =  total  head  manometer  difference  in  feet 

Ahg  =  static  head  manometer  difference  in  feet 

Az  =  difference  between  total  and  static  pressure- 
side  manometer  readings,  in  feet 

Because  of  small  variations  in  values  of  U0  for  dif- 
ferent sets  of  data,  it  was  decided  to  normalize  u  with  respect 
to  the  value  of  u  at  y/t  =  +1.5  in  every  case. 
C.  Turbulence  Intensity,  /y/u»^ 

Following  the  development  set  forth  in  Section  IV,  the 
turbulence  intensity, y    u'2  ,  was  calculated  from  the  rms 
turbulence  voltage  measurements  using  the  expressions 
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-  V, 


V   ui  2  =  JSS 
u 


(29) 


GaGpEtg 

Turbulence  Intensity  was  normalized  with  respect  to  the 
value  of  u  obtaining  at  the  point  of  observation. 

D.  Static  Pressure,  p. 

Static  pressure  was  calculated  from  the  static  head 
tube  manometer  reading  for  each  measurement  location.,   It 
was  normalized  in  terms  of  a  pressure  coefficient  P  "  Po 

where  p0  is  the  free  stream  pressure  and  U0  is  the  free 
stream  velocity. 

E.  Normalized  Turbulence  Energy  Density 9  E(n)/  u'2. 

Normalized  turbulence  energy  distribution  density  was 
calculated  directly  from  the  thermocouple  voltmeter  reading 
and  the  reading  of  the  vibration  analyzer  at  each  observation 
frequency  from  the  expression: 


Mai  = L 


u 


•  2 


V"rms  0V 


Rd2(n) 
n 


(30) 


This  normalization  process  eliminates  variations  in 
absolute  magnitude  of  energy  density  due  to  variations  in 
u  existing  over  the  entire  investigation. 
F,  Profile  Drag  Coefficient,  cDp. 

The  profile  drag  coefficient  was  calculated  at  longi- 
tudinal  stations  x/t  =3.25  and  x/t  =  9 .75*  using  the 
expressions 


°Dp  = 


Sa   [i-a-i- 


Po 


u 


u, 


|PU02 


(30 


'o      wo 
where  the  integration  was  performed  numerically  by  the 

trapezoidal  rule  between  the  limits  y/t  =  -1.5  and 

y/t  =  +1.5. 
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Applicability  of  this  expression  to  the  early  wake  was 

checked  using  the  Jones  Method1-  JJ  where: 


%  = 2 


j 


^(t-ywy&     (32) 


where:  g2  =  total  pressure  at  the  measuring  station 
P2  =  static  pressure  at  the  measuring  station 
p  =  static  free  stream  pressure 
q  =  total  pressure  minus  static  pressure  in 
the  free  stream. 
Integration  was  performed  numerically  as  above. 

It  should  be  noted  that  the  drag  coefficient  in  terms 
of  wetted  surface  is  given  by: 
CDp  =  16  C^ 

for  a  chord-to-thickness  ratio  of  16  as  is  the  case  for 
these  test  plates'-  ^  . 
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VII.  RESULTS  OF  THE  INVESTIGATION 

A.   Wake  Structure. 

The  results  of  the  study  of  the  wake  structure  of  the 
three  test  plates  are  presented  graphically  as  Appendix  B 
to  this  report.   The  salient  points  evidenced  by  these  plots 
will  be  discussed  in  this  section. 

Figs.  B-l  through  B-3  are  summary  plots,  showing  for 
each  of  the  three  test  plates  all  of  the  quantities  investi- 
gated in  the  wake  surveys.   Necessarily  drawn  to  a  reduced 
scale,  these  plots  are  intended  merely  to  indicate  in  one 
location  the  general  pattern  of  the  normalized  mean  velocity 
profiles  and  longitudinal  distributions,  the  transverse  and 
longitudinal  distributions  of  normalized  turbulence  intensity, 
and  the  transverse  and  longitudinal  distributions  of  normal- 
ized static  pressure. 

Figs.  B-ij.  through  B-6  show  for  each  individual  plate 
the  development  of  the  dimensionless  velocity  profile. 
Values  of  u  are  normalized  in  each  case  with  the  faired 
value  of  u  obtaining  at  station  y/t  =  +1.5.   Profiles  are 
shown  for  positive  values  of  y/t  only.   (This  is  the  case  for 
all  remaining  wake  plots  with  the  exception  of  the  wake 
symmetry  plots.)  None  of  the  plates  show  any  distinct  in- 
dication of  rapid  wake  spreading.   If  the  limits  of  the  wake 
are  arbitrarily  supposed  to  be  the  value  of  y/t  for  which 
^>!"ffi<  =  0,  it  is  seen  that  for  Plate  A  the  limits  of  the 
wake  gradually  increase  from  about  y/t  =  1.0  to  y/t  =  1.2 
between  the  four  longitudinal  stations.   For  Plate  B  the  wake 
spreads  from  about  y/t  =  0.8  to  y/t  =  1.2$  while  for  Plate  C 
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the  wake  spreads  from  about  y/t  =  1.0  to  y/t  =  I.I4..   The 
velocity  defects  for  Plates  A  and  C  at  x/t  =  3.25  are  seen 
to  be  markedly  greater  than  for  Plate  C,  and  the  decay  of  the 
velocity  defect  Is  more  rapid  for  these  two  plates.  The 
transverse  distributions  of  velocity  defect  are  essentially 
Gaussian. 

The  Increased  velocity  defect  at  the  early  stations  for 
Plates  A  and  C  Is  illustrated  further  in  Figs.  B~7  through 
B-10,  which  present  a  comparison  of  the  dimensionless  velocity 
profiles  at  each  of  the  longitudinal  measuring  stations. 
Close  to  the  trailing  edge  both  of  these  plates,  which  have 
relatively  sharp  separation  points,  have  large  defects,  al- 
though at  x/t  =  9.75  all  three  profiles  are  essentially  the 
same.   Plate  A,  with  its  90°  reentrant  angle  trailing  edge, 
is  seen  to  be  Intermediate  between  the  other  two  plates. 

Figs.  B-ll  through  B-ll+  present  the  transverse  dis- 
tribution of  the  normalized  turbulence  intensity  for  all 
three  plates  at  each  of  the  longitudinal  stations.  The  wake 
behind  Plate  C  is  seen  to  be  the  most  turbulent  in  the 
centerline  areaj  that  of  Plate  B  is  the  least  turbulent. 
Again  Plate  A  is  intermediate.  Turbulence  intensity  decays 
rapidly  downstream  and  at  x/t  =  9«7f>  the  intensity  levels 
are  essentially  equal  and  approaching  a  constant  value  across 
the  wake,  indicating  a  decay  to  isotropic  turbulence.  Tur= 
bulence  intensities  appear  to  approach  the  free  stream  level 
by  y/t  =  1.5  at  each  station  for  each  plate,  and  the  free 
stream  level  of  turbulence  intensity  is  seen  gradually  co 
increase  as  energy  is  diffused  from  the  wake  into  the  stream. 
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Fig.  B-15  through  B-18  illustrate  the  transverse  dis- 
tributions of  static  pressure  coefficient  in  the  wake  of  the 
three  plates  at  each  of  the  longitudinal  stations.   Plates  A 
and  C  with  their  more  sharply  defined  separation  points  are 
seen  to  have  extensive  underpressure  regions  at  x/t  =  1.5* 
Plate  B  with  its  rounded  trailing  edge  does  not.   As  we  move 
downstream  the  Plate  B  distribution  becomes  partially 
negative,  and  all  three  plots  move  together.   At  x/t  =9.75 
(see  Fig.  B-18)  the  distribution  plot  for  Plate  A  is  seen  to 
become  markedly  negative  again.  This  particular  presentation 
is  not  substantiated  by  the  investigation  of  longitudinal 
pressure  coefficient  distribution  illustrated  in  Figs.  B-21 
and  B-22,  and  may  be  due  to  data  error  or  computation  error. 
However,  the  possibility  that  this  plot  is  correct  must  not 
be  entirely  discounted  in  view  of  the  unusual  trailing  edge 
configuration  of  this  plate.  The  Plate  A  plot  in  Fig.  B-l£ 
shows  sharp  excursions  which  may  be  not  entirely  due  to 
scatter.  That  the  sharp  drop  in  pressure  coefficient  at 
y/t  =  0.5  may  not  be  due  to  data  scatter  is  suggested  by  the 
fact  that  Plate  B  also  indicates  a  smooth  decrease  in  this 
region.   It  will  be  seen  later  that  Plate  A  exhibits  a  very 
sharp  longitudinal  pressure  coefficient  gradient  in  this 
region  and  that  reverse  flow  is  also  observed. 

Fig.  B-19  presents  longitudinal  distributions  of  normal- 
ized mean  local  velocity  at  the  centerline  and  at  y/t  =  l.fj. 
In  these  plots  u  is  normalized  with  respect  to  free  stream 
velocity,  U0.  Reverse  flow  was  observed  for  Plate  A  in  the 
region  from  x/t  ■  0.75  to  about  x/t  =  1.25»  and  a  sharp 
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longitudinal  velocity  gradient  exists .   Reverse  flow  was  also 
observed  for  Plate  C  close  to  the  trailing  edge. 

Fig.  B-20  shows  longitudinal  distributions  of  turbulence 
intensity  for  the  three  plates  at  the  wake  centerline  and  at 
y/t  =  1.5.  Again  Plate  C  indicates  the  highest  level  of  wake 
turbulence  and  Plate  B  the  lowest °,    Plate  A  is  intermediate. 

FigSo  B-21  and  B-22  illustrate  the  longitudinal  static 
pressure  distributions  for  the  wake  centerline  and  at  y/t  = 
l.£  respectively.  Fig.  B-21  clearly  shows  the  large  under- 
pressure in  the  wake  close  to  the  trailing  edge  of  Plate  A. 
Plate  C  also  demonstrates  considerable  early  underpressure, 
while  Plate  B  underpressure  is  relatively  minor  and  very 
rapidly  becomes  an  overpressure  for  a  short  period,  as  does 
the  plot  for  Plate  A„  The  centerline  pressure  for  Plate  C 
is  seen  to  remain  negative  throughout  the  longitudinal  range. 

Figs.  B-23  through  B-25  Illustrate  the  essential 
symmetry  of  the  wakes  behind  the  three  plates.  Although  some 
data  scatter  is  seen,  it  appears  that  the  exact  centerline  of 
the  flow  was  slightly  on  the  positive  side  of  the  established 
axis,  fortuitously  insuring  that  the  flow  centerline  was 
actually  observed  during  the  half-wake  surveys  as  well  as  in 
the  full-wake  surveys. 
B.   Profile  Drag  Coefficients. 

Profile  drag  coefficients  were  calculated  for  the  three 
test  plates  at  longitudinal  stations  x/t  =  3.25  and  x/t  =  9°75' 
Calculations  were  initially  made  using  Eq.  (31),  and  then 
using  the  Jones  method  as  indicated  in  Eq.  (32).  Results  of 
these  calculations  are  listed  as  follows  in  Table  II. 
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Profile  Drag  Coefficients 
PLATE     x/t       °Dp  [Eq.  (31)]       GDp  [Eq.  (32)] 


A 

3.25 

0.443 

0.141 

A 

9.75 

0.378 

0.171 

B 

3.25 

0.519 

0.351 

B 

9.75 

0.318 

0.195 

0 

3.25 

0.528 

0.196 

0 

9.75 

0.393 

0.308 

The  above  values  of  °Dp  computed  from  Eq.  (31 )  are 
corrected  for  the  longitudinal  pressure  gradient  by  alge- 
braically subtracting  from  the  momentum  term  a  pressure 
coefficient  defect  term.   This  term  was  obtained  by  numer- 
ically integrating  the  defect  of  the  pressure  coefficient 
plot  where  the  defect  was  taken  as  zero  at  y/t  =  +1.5.   The 
results  compare  reasonably  with  the  value  of  °Dp  =  0.537 
predicted  by  Eq.  (16).   The  results  yielded  by  Eq.  (32)  are 
considered  incorrect. 
0,   Vortex  Shedding  Frequency. 

Vortex  shedding  frequencies  were  measured  by  observing 
the  peak  responses  in  the  normalized  energy  density  spectra 
for  the  three  test  plates  (see  Appendix  C).   These  shedding 
frequencies  were  observed  most  strongly  at  y/t  =  0.59  in  line 
with  the  edge  of  the  test  plate.   The  double  shedding  frequency 
is  more  clearly  noticeable  in  the  spectra  taken  at  y/t  =  0, 
on  the  centerline0   In  some  cases  the  triple  shedding  frequency 
is  also  noticeable. 

Plate  thickness  Strouhal  numbers  were  calculated  for 
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each  of  the  three  plates,  based  upon  averaged  flow  velocities 
and  observed  shedding  frequencies;  the  results  of  these  cal- 
culations are  presented  in  Table  III. 

In  order  to  compare  the  results  of  these  calculations, 
it  was  decided  to  normalize  the  Strouhal  numbers  with  respect 
to  Plate  C.   It  was  reasoned  that  for  a  small  Reynolds  number 
range  the  Strouhal  number  might  be  expected  to  remain 
essentially  a  constant.   Since  variations  in  Strouhal 
number  must  be  due  to  variations  in  the  transverse  separation 
between  vortices,  Plate  C  was  selected  as  a  basis  for  com- 
parison since  the  square  trailing  edge  provides  distinct 
separation  points  in  comparison  with  either  of  the  other  two 
test  plates.   It  was  then  assumed  for  the  sake  of  comparison 
that  the  transverse  separation  between  vortices  is  indeed 
equal  to  plate  thickness  for  Plate  C.   The  Plate  C  Strouhal 
number  was  then  taken  as  being  the  constant  value  for  all 
three  plates,  neglecting  small  variations  due  to  slightly 
different  plate  thickness  Reynolds  numbers,  and  this  Strouhal 
number  was  considered  a  vortex-separation  Strouhal  number, 
Sb  =  nb/U.  From  this,  the  relative  values  of  vortex  separ- 
ation were  calculated  for  Plates  A  and  B.   These  results  are 
also  listed  in  Table  III,  and  their  physical  significance 
is  illustrated  in  Fig.  VII. 
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TABLE 

III 

Strouhal  Numbers 

and  Relative 

Transverse  Vortex 

Sp 

aclng 

PLATE 

o.  _  nt 
St  -  -jf- 

0.237 

n 

U 

0  _  nb 
Sb  "  IT 

0.223 

b 

A 

109.87 

9.66 

0.236" 

B 

0.2^1 

llk.87 

9.53 

0.223 

0.222" 

C 

0.223 

103.00 

9.63 

0.223 

0.250" 
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FIGURE  VII.  Relative  Transverse  Vortex  Spacing 

The  qualitative  representation  of  vortex  spacing  in  Fig. 
VII  illustrates  expected  behavior  in  the  wake  of  Plate  B; 
the  separation  points  move  toward  the  centerline  of  the  plate 
for  a  semi-cylindrical  trailing  edge.  When  it  is  realized 
that  vortices  do  not  appear  at  the  separation  points  as  dis- 
crete entities,  but  are  gradually  formed  from  vorticity 
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which  diffuses  out  of  the  boundary  layer,  even  the  behavior 
of  the  vortex  locus  of  Plate  A  appears  plausible.   A  "cavity" 
of  relatively  large  under-pressure  exists  immediately  behind 
the  reentrant  wedge,  as  was  pointed  out  in  Section  VII, A.  and 
illustrated  in  Appendix  B,  and  the  flow  would  be  expected  to 
be  drawn  toward  the  centerline  of  the  wake.   The  loci  of 
vortices  are  illustrated  as  gradually  diverging  from  the 
wake  centerline  in  all  three  cases  in  conformance  with  the 
indication  of  wake  spreading  seen  in  Appendix  B. 

Observed  shedding  frequencies  were  compared  with  the 
frequencies  predicted  from  Eq.  (19).   This  comparison  is  pre- 
sented below  in  Table  IV.   Again,  exact  agreement  is  not 
realized,  but  a  reasonable  approximation  is  indicated. 

TABLE  IV 

Comparison  of  Actual  and 
Predicted  Shedding  Frequencies 


PIiATE 

vactual 

n,r 
vpredicted 

Eq.  (19) 

A 

110 

113.6 

B 

115 

112.0 

0 

103 

113.1 
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D.   Turbulence  Energy. 

The  results  of  the  investigation  of  the  turbulence 
energy  contained  at  various  stations  in  the  early  wake  behind 
the  three  test  plates  are  presented  graphically  in  Appendix  0. 

Since  minor  variations  existed  in  the  flow  conditions 
obtaining  for  the  individual  locations,  it  was  felt  that  no 
precise  comparison  could  be  made  of  absolute  magnitudes  of 
energy  density.   To  this  end,  the  spectra  were  normalized  as 
indicated  earlier.   The  data  presented  on  these  plots  is 
independent  of  the  local  velocity  at  the  point  of  measure- 
ment. 

The  area  under  these  spectrum  curves  is  dimensionless 
and  is  a  measure  of  the  total  energy  present  at  a  particular 
location  in  the  wake.   Integration  of  the  energy  density 
spectrum  from  zero  to  infinity  should  properly  result  in  a 
value  of  unity;  this  is  not  the  case  for  any  of  the  spectra 
obtained.   This  Indicates  that  a  certain  percentage  of  the 
total  energy  present  is  either  not  measured  or  not  accounted 
for  in  some  manner.   Possible  explanations  for  this  fact  are: 

(1)   The  limits  of  integration  are  the  finite  limits 
of  the  wave  analyzer.   Clearly  some  energy  must  be  present 
below  2.5  cps  and  above  750  cps.   Subsequent  measurement  pf 
one  of  the  spectra  with  a  General  Radio  Type  1 554— A  vibration 
analyzer  indicates  a  considerable  amount  of  energy  present 
at  about  1400  cps,  which  is  possibly  attributable  to  Helm- 
holtz  resonance  effects.  (See  Appendix  D. )  It  is  not  believed 
that  any  significant  amount  of  energy  is  present  below  the 
low  frequency  limit  of  the  wave  analyzer  because  of  the 
small  frequency  Increment  under  consideration. 
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FIGURE  VIII  Comparison  of  Wave  Analyzer  Energy  Density  Spectra 


50 

(2)   The  spectral  measurements  made  with  the  wave  anal- 
yzer are  normalized  with  respect  to  measurements  made  with  a 
true  rms  voltmeter.   It  is  possible  that  the  thermocouple 
voltmeter,  being  an  all-pass  energy  measuring  device  and 
Indicating  a  true  picture  of  the  total  energy  present  at  a 
given  location,  measures  components  which  the  wave  analyzer 
meter  does  not.   The  wave  analyzer  is  a  tunable  voltmeter 
designed  to  measure  Fourier  components  of  vibration,  and  some 
energy  components  may  not  be  recognized  by  it.   During  the 
course  of  the  energy  spectrum  measurements,  simultaneous 
readings  were  taken  with  the  Ballantine  VTVM.   On  the  assumption 
that  the  foregoing  might  be  correct,  the  energy  spectra  were 
also  normalized  with  respect  to  the  Ballantine  voltage,  since 
it  was  reasoned  that  this  instrument  is  similar  to  the  meter 
installed  in  the  wave  analyzer.  The  results  of  these  calculations 
which  are  not  plotted,  still  Indicate  an  energy  defect,  but  a 
much  smaller  and  more  reasonable  one. 

Further  investigation  will  be  necessary  to  explain  this 
discrepancy  fully  before  any  absolute  magnitudes  can  be 
assigned  to  energy  data. 

Because  of  the  foregoing,  no  magnitude  comparison  of 
discrete  energy  represented  by  the  vortex  shedding  frequencies 
could  be  made.   In  order  to  examine  the  vortex  energy  as  a 
function  of  trailing  edge  configuration  then,  the  ratio 
Ap/Ag  was  studied,  where  A  is  the  area  under  a  discrete 
shedding  frequency  Cor  double  shedding  frequency)  peak  and 
Ag  is  the  total  area  under  the  spectral  plot  obtained.   This 
ratio  is  then  a  measure  of  the  percentage  of  energy  present 
in  the  vortex  shedding  phenomenon  for  each  of  the  spectra. 
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The  ratio  Ap/Ag  is  greatest  at  any  given  location  for 
Plate  B  and  least  for  Plate  0;  Plate  A  is  again  intermediate. 
This  fact  was  noted  for  both  the  fundamental  and  double 
shedding  frequencies.   The  ratio  Ap/As  shows  a  linear  decay- 
downstream  for  both  Plates  B  and  0.   The  value  for  Plate  A 
increases  to  a  maximum  in  the  vicinity  of  x/t  =  3.25  and 
then  exhibits  essentially  a  linear  decay  with  x/t. 

Two  independent  investigations  of  the  ambient  level 
of  turbulence  in  the  empty  test  section  (with  no  test  plate 
installed)  were  conducted.   The  normalized  energy  density 
spectra  for  these  Investigations  are  presented  in  Pig.  0-15. 
Both  spectra  were  taken  at  the  tunnel  location  corresponding 
to  x/t  =  3.25  and  y/t  =  0.   For  each  investigation  the  value 
of  u  measured  by  the  total  and  static  head  tubes  was  found  to 
be  9.52  fps.   Turbulence  intensities  for  these  measurements 
were  0.0352  and  0.0366,  considerably  lower  than  the  value 
of  0.068  found  by  Huval  when  operating  at  u  =  10o2  fps. 

Also  presented  in  Pig.  0-15  are  the  theoretical 
functions  suggested  respectively  by  Dryden  and  von  Karman: 


0(n)  =    A1        and  0(n)  = 


A« 


1  +  B^n2  (1  +  B22n2)5/6 

While  reasonable  agreement  is  obtained  with  either  function 
it  appears  that  the  experimental  data  are  better  represented 
by  the  Dryden  expression  in  this  case.   The  individual  plate 
wake  spectra  are  also  well  represented  by  this  expression  in 
the  high  frequency  region.  This  is  an  indication  of  an 
approach  to  isotropic  turbulence  in  this  region,  while  anisotropy 
is  indicated  in  the  low  frequency  plate  spectra  regions. 
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Pigs.  0-1  through  C-3  illustrate  the  downstream  develop- 
ment of  the  turbulence  spectra  for  the  individual  plates 
along  the  centerllne.   Figs.  C-4  through  0-6  present  the 
longitudinal  development  for  the  individual  plate  spectra 
in  line  with  the  edge  of  the  plate.   Pigs.  C-7  through  0-10 
compare  the  centerllne  spectra  for  the  three  test  plates 
at  each  of  the  four  longitudinal  stations,  and  Pigs.  C-11 
through  0-14  make  a  similar  comparison  for  the  spectra 
obtained  at  y/t  =  0.5.   Inspection  of  these  Figures  shows 
the  following  trends  which  are  worthy  of  note : 

Between  about  400  cps  and  750  cps  the  slope  of  all 
the  spectra  is  essentially  -2  on  the  log-log  plot.   At 
frequencies  between  the  double  vortex  shedding  frequency 
and  about  400  cps  the  slope  is  well  represented  by  -5/3. 

Discrete  energy  density  spikes  associated  with  vortex 
shedding  are  very  pronounced  at  x/t  =1.50  and  at  x/t  = 
3.25  for  both  y/t  =  0  and  y/t  =  0.5.   Spikes  associated 
with  the  double  shedding  frequency  are  pronounced  only  at 
the  centerllne. 

Considerable  scatter  exists  in  the  spectra  in  the  frequency 
range  from  2.5  cps  to  about  20  cps. 

The  discrete  frequency  energy  appears  to  diffuse  initially 
to  the  random  components  immediately  adjacent  in  frequency. 
These  components  in  turn  transfer  part  of  their  energy  to 
their  neighboring  components.   Thus,  as  the  flow  moves  down- 
stream the  discrete  energy  spikes  are  rapidly  reduced  by 
diffusion  of  their  energy  to  the  random  components  in  a  sort 
of  "relaxation"  process. 

As  the  flow  moves  downstream  from  x/t  =  1.50,  the 
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turbulence  spectra  indicate  transition  from  very  pronounced 
anisotropic  turbulence  superposed  upon  strong  discrete 
components  of  vorticity  toward  a  more  isotropic  condition 
of  statistical  equilibrium.   Apparently  the  high  frequency 
components  achieve  equilibrium  early  in  the  wake,  since  the 
portions  of  the  spectra  beyond  300  cps  do  not  appear  to 
undergo  much  change  as  the  flow  moves  downstream.   Hinze 
suggests  a  plausible  reason  for  this  early  achievement  of 
stability  in  the  high  frequency  end  of  the  spectrum.   He 
points  out  that  the  flow  within  the  smallest  eddies  is  no 
longer  turbulent,  but  viscous,  and  that  the  lower  limit  on 
turbulent  eddy  size  is  determined  by  dominance  of  molecular 
effects    '.   In  general,  the  severity  of  velocity  grad- 
ients increases  with  decreasing  eddy  size.   The  rate  of 
viscous  dissipation  of  energy  associated  with  large  velocity 
gradients  then  sets  the  lower  limit  on  eddy  size  in  accordance 
with  the  principles  of  conservation  of  energy. 

Numerical  integration  of  the  first  four  energy  density 
spectra  taken  for  Plate  0  resulted  in  very  large  energy  defects 
compared  to  all  the  other  spectra  taken.  Applying  the 
correction  suggested  above  as  a  possible  means  of  adjusting 
for  differences  in  measuring  instruments  (i.e.,  normalizing 
with  respect  to  the  Ballantine  voltage  recorded  instead  of 
the  voltage  reading  of  the  thermocouple  voltmeter)  only 
slightly  improved  the  situation.  For  these  spectra  Ag  still 
remained  less  than  50$  . 

Possible  explanations  for  this  large  discrepancy  are: 
(1)   Incorrect  wave  analyzer  calibration  constant.  This 
is  discussed  further  in  Appendix  A. 


54 

(2)  Abnormally  large  amounts  of  energy  associated  with 
frequencies  above  750  cps.   As  shown  in  Pig.  VIII,  an  extended 
spectrum  analysis  taken  with  the  Type  1 554— A  analyzer  indicates 
about  }0%   of  the  spectral  energy  above  750  cps.   However,  the 
energy  above  750  cps  would  have  to  be  considerably  larger  than 
this  to  account  for  the  discrepancy. 

(3)  Presence  of  an  air  bubble  in  the  probe  tip.   If 
the  natural  response  of  the  gage-water-column-probe-tip 
system  is  located  at  a  frequency  just  above  the  range  of  the 
wave  analyzer,  a  large  amount  of  energy  will  exist  at  this 
point  and  will  not  be  measured  by  the  analyzer.   It  will,  how- 
ever, be  measured  by  the  rms  voltmeter.   This  is  considered 
unlikely,  however,  since  two  separate  probe  tip  filling 
operations  were  Involved. 

Pig.  VIII  indicates  that  the  relative  magnitude  of  the 
four  Plate  C  spectra  in  question  are  probably  correct,  since 
the  pattern  is  essentially  duplicated  by  the  spectrum  obtained 
with  the  Type  1 55^-A  analyzer. 

In  order  to  determine  whether  the  probe  tip  itself  was 
subject  to  vibrations  in  the  highly  fluctuating  flow  field , 
the  probe  tip  opening  was  plugged  with  candle  wax.   The 
turbulence  probe  was  then  located  at  x/t  =3.25  and  y/t  =0.5 
in  the  wake  of  Plate  0.   With  a  flow  of  U0  =  9.55  fps,  a 
turbulence  intensity  of  0.0019  was  obtained.   With  the  KinTel 
amplifier  gain  increased  to  100,  the  entire  spectrum  was  scanned; 
response  was  negligible  at  all  frequencies.   With  the  flow  secured 
a  voltage  was  still  measured  with  the  thermocouple  voltmeter; 
this  background  signal,  referred  to  a  flow  rate  of  9.55  fps, 
corresponded  to  an  intensity  of  turbulence  of  0.0006. 
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A  graphical  presentation  of  the  decay  of  discrete  vortex 
energy  is  given  below  in  Fig.  IX: 
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FIGURE  IX  Decay  of  Vortex  Energy  Downstream. 
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VIII.   CONCLUSIONS 

A  study  of  the  graphical  presentation  of  results  included 

as  Appendices  B  and  C  indicates  the  following  trends  for  these 

4 
test  plates  at  a  plate-thickness  Reynolds  number  of  1.9  x  10  : 

(1)  Sharp  geometry  transitions  (well  defined  flow 
separation  points)  in  the  trailing  edge  create  a  higher  level 
of  turbulence  intensity  in  the  early  wake  than  do  gradual 
transitions.   Turbulence  intensity  diffuses  rapidly  out  of 
the  wake  at  a  rate  proportional  to  the  initial  magnitude  of 
the  Intensity  and  gradually  approaches  a  decay  law  propor- 
tional approximately  to  (x/t)~2. 

(2)  Strong  wake  underpressures  exist  in  the  early  wake 
behind  sharp  geometry  transition  trailing  edges. 

(3)  The  energy  density  spectrum  in  the  early  wake 

is  anisotropic,  particularly  close  to  the  trailing  edge  where 
longitudinal  and  transverse  pressure  gradients  are  greatest, 
and  particularly  for  low  frequencies. 

(4)  The  high  frequency  end  of  the  energy  density 
spectrum  in  the  early  wake  is  seen  to  approach  isotropic 
conditions.   This  is  also  true  of  the  ambient  tunnel  spectrum. 
E(n)  is  seen  to  vary  essentially  as  n"2  in  this  region. 

(5)  The  discrete  energy  represented  by  the  vortex 
shedding  frequency  is  seen  to  decay  downstream  in  a  manner 
similar  to  the  decay  of  turbulence  intensity.   The  rounded 
trailing  edge  geometry  favors  the  formation  of  vortices  of 
higher  energy  content. 

(6)  The  magnitudes  of  the  turbulence  intensity  found 
in  the  early  wake  close  to  the  trailing  edge  are  in  excess 
of  the  order  of  magnitude  originally  assumed  in  the 
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development  of  the  expression  for  turbulence  gage  response. 
These  magnitudes  must,  therefore ,  be  viewed  in  a  qualitative 
sense  only. 

(7)  A  significant  percentage  (as  much  as  50$)  of  the 
total  kinetic  energy  of  a  mixed  periodic-turbulent  wake  is 
associated  with  discrete  vortex  shedding  frequencies,, 

(8)  Energy  associated  with  vortex  shedding  frequencies 
increases  initially  in  the  case  of  the  plate  with  a  90° 

reentrant  angle  trailing  edge  (Plate  A),  and  then  decays 

(13) 

approximately  linearly  downstream.   Huval     found  similar 

behavior  In  the  early  wake  behind  a  60°  included  angle  trail- 
ing edge  geometry. 

(9)  Taking  the  amount  of  energy  associated  with  the 
discrete  shedding  frequencies  in  the  Immediate  vicinity  of 
the  trailing  edge  as  a  measure  of  the  energy  available  to 
excite  plate  vibrations,  extrapolation  of  the  trends  indicated 
in  Pig.  IX  Indicate  that  Plate  B  would  tend  to  vibrate  most 

and  Plate  A  would  vibrate  least.   This  is  in  agreement  with 

(14) 
results  reported  by  Ippen  et  al 

While  centerline  vortex  energy  close  to  the  trailing  edge 
is  greatest  for  Plate  A,  wake  centerline  forces  exerted  on  the 
trailing  edge  will  act  equally  on  either  side  of  the  reentrant 
angle;  response  to  such  forces  would  thus  inherently  be  damped, 
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IX.   RECOMMENDATIONS 

The  following  recommendations  are  made  for  the  conduct 
of  further  Investigation  of  the  early  wake  behind  flat  plate 
structures i 

(1)  The  predictability  of  vortex  shedding  frequency  and 
the  behavior  of  plate-thickness  and  vortex-spacing  Strouhal 
number  should  be  investigated  as  a  function  of  Reynolds 
number.   Vortex  spacing  should  be  investigated  by  traversing 
the  early  wake  with  the  turbulence  transducer  and  obtaining 
transverse  distributions  of  discrete  spectral  energy. 

(2)  The  behavior  of  profile  drag  coefficient  with  var- 
iations in  trailing  edge  geometry  should  be  the  subject  of 

a  separate  investigation.   Techniques  and  Instrumentation 
in  this  investigation  were  Inadequate  for  any  conclusive 
understanding  of  such  behavior. 

(3)  Particular  care  should  be  taken  to  duplicate  both 
flow  velocity  and  tunnel  test  section  pressures  for  absolute 
comparison  of  successive  runs. 

(4)  Tunnel  wall  spindle  cavities  should  be  plugged  with 
plates  drilled  to  accept  plate  mounting  spindles  in  order  to 
eliminate  the  possibilities  of  Helmholtz  or  stub  tube  reson- 
ance effects.  (See  Appendix  D„ )   Any  other  cavities  in  the 
tunnel  walls  in  the  vicinity  of  the  test  plate  or  transducer 
should  be  similarly  plugged  or  faired. 

(5)  Frequent  periodic  calibration  of  measurement  equip- 
ment and  components p  particularly  when  they  are  battery 
powered,  is  mandatory  when  accurate  data  is  to  be  obtained 
over  a  period  of  time. 
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(6)  Accurate  analysis  of  energy  density  spectra  requires 

accurate  calibration  of  analysis  equipment.   True  rms  indication 
of  both  total  energy  level  and  energy  density  as  a  function 
of  frequency  is  desirable.   The  same  instrument  or  accurately 
cross-calibrated  Instruments  should  be  used  to  measure  these 
quantities . 

(7)  Consideration  should  be  given  to  utilizing  frequency 
modulated  tape  recording  equipment  to  permit  determination 

of  autocorrelation  functions  for  the  energy  density  spectra. 
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APPENDIX  A 


CALIBRATION  PROCEDURE 
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.  1.   Turbulence  Pressure  Transducer. 

a.  General. 

The  turbulence  gage  used  In  this  Investigation  was  the 

fio] 

same  as  that  employed  by  Huval1    .   A  quasi-static  calibration 
procedure  similar  to  that  used  by  Huval  was  used.   However, 
the  procedure  was  simplified  as  a  result  of  repairs  to  the 
transducer  made  by  Hydrodynamics  Laboratory  personnel  which 
resulted  in  a  greatly  increased  gage  voltage  decay  time 
constant. 

Huval  found  the  turbulence  gage  sensitivity  to  be  63 
millivolts  per  psi,  a  value  substantially  lower  than  the 
theoretically  predicted  sensitivity  of  100  millivolts  per 
psi.  Furthermore,  his  calibration1-  ~>i   showed  the  gage  time 
constant  to  be  in  the  vicinity  of  0.6  second,  considerably 
lower  than  that  to  be  expected.   Subsequent  to  Huval 's  in- 
vestigation the  sensitivity  of  the  gage  decreased  markedly. 
Hydrodynamics  Laboratory  personnel  then  disassembled  the 
gage  and  found  a  small  amount  of  foreign  matter  on  the 
ceramic  face.   The  ceramic  was  carefully  cleaned  and  re- 
assembled.  Tests  performed  prior  to  the  present  investi- 
gation disclosed  that  the  sensitivity  of  the  gage  was  in- 
creased and  the  gage  time  constant  was  greatly  improved. 

b.  Description  of  Calibration  Procedure. 

Fig.  A-l  shows  schematically  the  turbulence  gage  cali- 
bration equipment  used  in  this  investigation. 
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FIGURE  A-1   Gage  Calibration  Schematic 

The  turbulence  gage  with  probe  tip  installed  was  in- 
serted in  the  tube  connected  to  the  side  of  the  cylindrical 
pressure  tank  as  shown.   A  differential  pressure  was  built 
up  in  the  mercury  manometer  by  means  of  the  hand-operated 
pressure  bulb.  The  clamp  in  the  manometer  pressure  tube  was 
set  to  give  sufficient  damping  constriction  to  reduce  pressure 
oscillations.   The  solenoid  valve  at  the  top  of  the  pressure 
chamber  was  opened  by  means  of  a  trigger  circuit  which  simul~ 
teneously  triggered  the  horizontal  sweep  on  the  oscilloscope. 
An  oscilloscope  camera  was  used  to  record  the  indicated  gage 
response  voltage.   See  Fig.  A-2. 

A  simpler  method  of  calibration,  eliminating  the  use  ©f 
the  oscilloscope  camera,  was  employed  for  the  majority  of  the 
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turbulence  gage  sensitivity  calibrations  in  this  investigation. 
In  this  alternative  method,  the  oscilloscope  sweep  rate  was 
increased  and  the  horizontal  sweep  selector  set  to  "recurring 
sweep."  The  vertical  deflection  of  the  horizontal  sweep  at 
time  of  pressure  release  was  measured  directly  on  the  face  of 
the  oscilloscope,  manually  recording  on  a  strip  of  paper. 
One  series  of  photographic  calibration  runs  was  made  immediate- 
ly after  a  series  of  runs  using  the  simplified  procedure. 
Reproducibility  between  the  two  methods  was  excellent,  both 
yielding  the  same  value  for  gage  sensitivity. 

c.  Development  of  the  Expression  for  Turbulence  Gage 


Sensitivity  and  Time  constants » 
The  following  definitions  will  be  used  in  this  develop- 
ment: 

K^g  =  turbulence  gage  sensitivity 

Kfcp  =  transducer-preamplifier  system  sensitivity 

Gp  =  preamplifier  gain 

Tv  =  voltage  time  constant 

Tp  =  pressure  time  constant 

Cv  =  1/TV 

Cp  =  1/Tp 

p   =  instantaneous  pressure  at  transducer 

pQ  =  initial  calibration  pressure 

V   =  instantaneous  voltage  at  oscilloscope 
As  indicated  by  Huval1  JJ,  the  sudden  release  of 
pressure  from  the  tank,  assuming  a  suitably  damped  system, 
may  be  represented  by  the  following  expressions: 
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p  =  p0e"CPt  .  for  t  >0 

P  =  P0  ,  for  t  <0  (A-l) 

The  incremental  change  in  voltage  generated  by  the  tur- 
bulence gage  is  the  sum  of  two  components: 

AV"i  =  -Kj-pAp    (pressure-generated  component) 

AV2  =  -CvVAt    (voltage  decay  component) 
Hence:     AV  =  AVX  +  AV2  =  -KtpAp  -  CvVAt  (A-2) 

Rewriting  Eq.  (A-2)  as: 

and  taking  the  limit  as  At  -*•  0,  we  obtain: 

Substituting  p  =  p0e~  P  into  Eq.  (A-3)  we  obtain 

H  =  KtpPoCpe"Cpt  "  CvV  (A^> 

Defining  VQ  =  K+-DP0»  we  see  that: 

H  =  Vpe"Cpt  -°vV  (A-5) 

Multiplying  both  sides  of  Eq.  (A-5)  by  e°vt: 

eCvt  dV  +  CyVeCvt  =  VoCp6 ( Cv-Cp ) t  ( A_6  j 


Then: 


&   (Ve^t)  =  VoCpe-V  (A-7) 


Integrating: 


VeCvt  .  VoCp  e«V.Cp)t  +  (A.6) 
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where  C2  is  the  constant  of  integration.   From  the  initial 
condition  that  V  =  0  at  t  =  0,  we  find  that: 

C2  =  -  ^£e_  (A-9) 

Then: 

V  =  M>_  fe-Cp*  -e-cvt)  (A.10) 

In  analysis  of  the  oscilloscope  calibration  photo- 
graphs, it  is  convenient  to  express  V  in  terms  of  Vra,  the 
maximum  voltage  occurring  at  time  tm.   (See  Fig.  A-2). 

From: 


dV 

vdT 


*  =  tm  =  ° 
m 


we  obtain: 


e-Cptm  =  %L  e"0**"  (A-ll) 

°P 


Also: 


Vm  =  ^°CP   (  e-%  -e-^m)  (A-12) 

Cv-Cp  \  / 

Substituting  Eq.  (A-ll)  into  Eq.  (A-12)  and  simplifying: 

Vm  =  T..-°v*-  =  KtpP0e-W*T  (A13) 

Solving  for  K^p  in  terms  of  Vm,  we  obtain: 

Ktp  -  —  e 

Before  going  on  to  the  development  of  expressions  for 
the  time  constants  Tv  and  Tp,  it  should  be  noted  from  the 
foregoing  expression  that  if  tm  <<^  Tv,  then: 


Ktp*$* 


P0  (A-15) 


This  yields  a  means  of  obtaining  an  adequate  approx- 
imation for  K.   without  accurate  knowledge  of  either  T 
or  Tp. 

Substituting  Eq.  (A-13)  into  Eq.  (A-10),  the  following 
expression  for  V  may  be  found: 
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v  =  T^-   .-<^> --%-£> 


rm 


(A-16) 


If  it  can  be  assumed  that  Tp  «  Tv  and  tm«Tv,  the 

second  term  within  brackets  in  Eq.  (A-16)  becomes  negligible 

with  respect  to  the  first  term  for  t»  tm. 

Then:  .  . 

-  ^  m 

(A-17) 


V   Vme 


;  for  t»  t 


m 


Hence: 


t-t 


▼  =  ln(Vnffy)  ->   evaluated  at  t»tm        (A-18) 


To  obtain  a  means  of  evaluating  Tp,  evaluating  dV/dt 
from  Eq.  (A-16)  above  at  t  =  0  results  in: 


dV 

dT 


t  =  0 


_  vm 

-  TtT- 


and: 


Tp  "  [dv/dtj 


;   evaluated  at  t  =  0  (A-19) 


This  allows  evaluation  of  Tp  from  knowledge  of  the 
initial  slope  of  the  V  versus  t  plot. 

d.  Results  of  Pressure  Transducer  Calibration. 


A  series  of  seven  photographic  calibration  runs  were 
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made  using  various  values  of  initial  pressure  pQ  and  using  an 
oscilloscope  horizontal  sweep  rate,  calibrated,  of  5>00  milli- 
seconds per  inch.   A  DuMont  Oscillograph-Record  camera  Type 
297 »  using  an  f-stop  of  f/5>.6  with  a  5-second  exposure,  the 
sweep  being  triggered  2  seconds  after  the  shutter  was  opened. 
Polapan  3000  Type  Ij-7  film  was  used.  The  time  to  reach  maximum 
voltage,  t  ,  was  about  30  milliseconds  on  each  run,  as  com- 
pared with  a  value  of  about  $0   milliseconds  reported  by  Huval. 

(1)  Pressure  Decay  Time  Constant,  Tp. 

Eq.  (A-19)  was  used  to  determine  Tp.  Values  ranged 
from  11  to  18  milliseconds,  with  five  of  the  seven  values 
being  in  the  range  of  12.5  to  11+.  milliseconds.  An  average 
value  of  Tp  =  13.7  milliseconds  was  obtained, 

(2)  Voltage  Decay  Time  Constant,  Tv. 

A  plot  of  ln(Vm/V)  versus  (t-tm)  was  used  to 
evaluate  Tp.  See  Eq„  (A-18).  A  range  of  values  from  2*1+ 
seconds  to  11. I4.  seconds  was  obtained,  with  five  of  the  seven 
runs  in  the  range  l+.lj.  to  9«0  seconds. 

Mubh  of  this  scatter  was  due  to  inability  to  obtain  an 
absolutely  stable  reference  voltage  on  the  oscilloscope  prior 
to  triggering  the  pressure  release.  Also,  extreme  sensitivity 
of  the  gage  and  preamplifier  to  small  unavoidable  disturbances 
during  the  2-second  period  of  the  horizontal  sweep  may  have 
contributed  to  the  scatter. 

The  two  extreme  values  were  disregarded,  resulting  in 
an  average  from  the  remaining  five  runs  of  Tp  =  6.1j.  seconds. 

However,  even  using  the  smallest  value  of  T  and  the 
largest  value  of  Tp,  it  is  seen  thats 
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m  +■ 

JB.  =  0.0075      and      wr  =  0.0125 

This  shows  that  the  approximation  of  Eq.  (A-l£)  is 
valid  for  determining  Ktp  and  Justifies  using  the  simplified 
procedure  described  in  Section  A.l.b  above. 
(3)  Turbulence  Gage  Sensitivity. 
Turbulence  gage  sensitivity  was  determined  from 
Eq.  (A-15)  and  the  relationship: 

Ktg  =  W% 

Seven  separate  series  of  calibration  runs  were  performed 
during  the  course  of  this  investigation  between  22  February 
and  11  April  1961.  The  dates  of  the  runs  are  presented  in 
the  legend  of  Pig.  A-3. 

During  the  observations  made  on  6  April,  the  KinTel 
amplifier  was  introduced  into  the  calibration  circuit  to 
provide  higher  output  voltage  levels.  Amplifier  gain  setting 
was  30.   It  was  felt  that  this  procedure  would  serve  as  a 
further  check  on  the  sensitivity  of  the  entire  measuring 
system.  Due  to  the  quasi-static  nature  of  the  calibration 
procedure,  the  DC  filter  between  the  preamplifier  and  the 
amplifier  was  bypassed.   The  average  value  of  K.   determined 
for  these  runs  was  about  20°/o  lower  than  the  February  cali- 
bration. 

As  a  check  on  possible  time-deterioration  of  sensitivity 
the  February  procedure  was  duplicated  on  11  Aprilj  these  runs 
produced  a  value  of  K.   about  \\°/o  lower  than  the  original 
February  value.   It  was  concluded  that  the  amplifier  muss 
have  attenuated  the  quasi-steady  gage  response. 
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FIGURE  A-3  Transducer-Preamplifier  Sensitivity  Calibration 
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The  amplifier  had  already  been  calibrated  using  an  AC 
signal  source  and  found  to  have  excellent  gain  response.   The 
amplifier  was  now  calibrated  using  a  DC  dry  cell  through  a 
potentiometer  as  a  signal  source.   With  gain  setting  of  30, 
actual  gain  was  measured  to  be  equal  to  about  27.   When  the 
ratio  of  gain  factors  was  applied  to  the  gage  calibration 
constants  obtained  with  the  amplifier  in  the  circuit,  the 
same  h°/o   variation  from  the  initial  calibration  sensitivity 
was  obtained.   It  is  considered  probable  that  this  variation 
is  due  to  data  scatter. 

Finally,  an  average  value  of  K.   =16.0  mv/ft.  of  water 
=  36.9  mv/psi  was  obtained.   This  corresponds  to  a  value  of 
K.   =  l+O  mv/ft.  of  water  =  92 •!+.  mv/psi,  as  compared  to  the 
theoretical  prediction  of  K.   =  100  mv/psi  reported  by 

r^2i 

Perkins  and  EaglesonUJ   . 
2.   Wave  Analyzer. 

Both  wave  analyzers  used  in  the  investigation  were 
calibrated  using  a  "pure  tone"  sinusoidal  input  signal  from 
the  Hewlett-Packard  audio  oscillator.   Rms  input  voltages 
were  measured  with  the  Ballantine  VTVM.   All  calibrations  were 
performed  with  the  vibration  analyzer  selectivity  in  the 
"sharp"  position  and  with  vibration  analyzer  sensitivity  set 
at  10.   (All  data  were  taken  with  these  settings  also.) 

a.  General . 

As  indicated  by  Perkins  and  Eagleson*-^   ,  the  General 
Radio  Type  762-B  Vibration  Analyzer  output  voltmeter  reading 
squared  ife  proportional  to  the  total  energy  level  of  the 
frequency  components  of  the  input  signal  within  the  narrow 
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band  to  which  the  analyzer  Is  tuned. 


b.  Wave  Analyzer  Response 


r,~   j. An 

E(n)dn  (A-20) 


_  An 
nC 


For  an  input  signal  with  a  continuous  spectrum  of 
frequencies,  if  E(n)  may  be  considered  constant  over  the 
narrow  bandwidth  of  the  wave  analyzer,  (A-20)  becomes? 

\t  =  E(ncMn  =  Ci  E(nc)nc  (A-21) 

From  the  statement  made  in  paragraph  (a)  above,  there 
results: 

Ebw  =  C2Rd2  (A-22) 

Therefore,  combining  (A-21)  and  (A-21)  °. 

E(nc)  =  t^-^  (A-23) 

°v   nc 

cl 
where  Cy  =  ^ 

c.  Wave  Analyzer  Calibration  Procedure. 

The  purpose  of  the  wave  analyzer  calibration  is  to  de- 
termine the  constant  C  in  (A-23)  • 

One  means  of  obtaining  this  constant  is  to  set  the  wave 
analyzer  at  a  given  frequency  n-,,  then  to  apply  an  input 
sinusoidal  signal  from  a  sine°wave  generator.  The  input 
signal  frequency  is  increased  in  small  Increments  from  a 
value  (well  below  n, )  where  negligible  wave  analyzer  response 
is  noted,  to  a  value  well  above  n-,  where  again  negligible 
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analyzer  response  results.  The  Input  rms  voltage  ei#  as 
measured  by  an  A.C.  rms  voltmeter,  should  be  kept  constant 
throughout  this  procedure  at  a  value  predetermined  to  give 
nearly  full  scale  deflection  at  maximum  wave  analyzer  response, 
e.  is  recorded  and  constantly  monitored.  At  each  input  fre- 
quency n,  n  and  Rd  are  recorded.  Then  Cv  may  be  determined 


from  the  relationships 

1 


Cv  = 


^-     dn  (A-2W 

nc 


The  foregoing  procedure  should  be  repeated  for  several 
values  of  n  over  the  range  of  frequencies  which  are  of 
interest.  By  this  means,  any  frequency  dependency  of  C  may 
be  detected. 

An  alternative  procedure  (introducing  a  small  error) 
is  to  set  the  input  frequency  at  a  constant  value  n.,,  then 
shift  the  wave  analyzer  frequency  in  small  increments  from 
the  point  of  negligible  response  on  the  lower  side  of  n,  to 
the  point  of  negligible  response  on  the  upper  side  of  n, » 
Then  Cy  may  be  approximated  from  the  relationships 


cv 7 

ei 


cod   ^ 

M_  dno  (A_25) 


The  latter  method  was  used  throughout  this  investigation, 
due  to  instability  of  the  sine  wave  generator  output  voltage 
introduced  by  shifting  frequencies.  Numerical  integration 
was  used  to  approximate  the  value  of  the  integral. 
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d.  Results  of  Wave  Analyzer  Calibrations. 

It  was  noted  during  this  investigation  that  the  wave 
analyzer  calibration  constant  was  dependent  upon  plate  vol- 
tage (supplied  by  two  \\$   volt  batteries). 

This  was  true  even  when  a  check  of  plate  battery  voltages 
showed  the  plate  voltage  to  be  within  the  satisfactory 
operating  range.  Therefore,  frequent  wave  analyzer  calibration 
is  necessary  in  order  to  obtain  accurate  results.   Particular- 
ly, the  wave  analyzer  must  be  recalibrated  whenever  the 
plate  supply  batteries  are  changed. 

It  is  possible  that  large  discrepancies  in  the  areas  of 
four  energy  spectra  for  Plate  C  may  be  attributed  to  failure 
to  calibrate  the  wave  analyzer  immediately  after  a  change  of 
plate  supply  batteries. 

On  March  28,  1961,  wave  analyzer  plate  supply  batteries 
were  replaced.  The  turbulence  spectrum  measurements  were 
commenced  and  all  eight  energy  spectra  for  Plate  C  were 
taken  between  March  28  and  March  31 »  when  the  wave  analyzer 
was  calibrated.  The  wave  analyzer  constant  was  found  to  be 
in  good  agreement  with  the  two  previous  values  determined. 
However,  it  was  later  discovered  that  the  integrals  of  the 
energy  density  spectra  were  very  low  for  the  four  spectra 
measured  on  March  28  and  March  29  compared  to  those  taken  on 
March  30,  and,  in  fact,  compared  to  all  the  spectra  taken  for 
the  other  two  plates.  At  present,  lacking  a  better  reason 
for  this  discrepancy,  it  is  conceivable  that  the  wave  analyzer 
constant  was  very  low  immediately  after  the  battery  replacement 
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for  some  unexplainable  reason*   In  evaluation  of  this  possi- 
bility, it  must  be  noted  that  a  very  low  value  compared  to 
any  previously  obtained  would  have  been  necessary  to  explain 
the  discrepancy,  and  the  constant  would  have  increased  with 
analyzer  operating  time,  contrary  to  the  trend  previously 
experienced. 

The  following  values  of  C  were  obtained  during  this 
investigation.  Only  average  values  are  listed,  since  each 
calibration  showed  little  frequency  dependence  of  the  constant 


PERIOD  Cv   (— ttt—  ) 
~~ —  volt^-cps 

3-6-61  to  3-10-61  5-12  x  105 
3=10-61  to  3-28-61  5.^1  x  lo£ 
3-28-61  to  3-31-61     5. 31  x  105 

3.  Amplifier  and  Preamplifier. 

The  electrometer  preamplifier  and  the  KinTel  amplifier 
were  both  calibrated  using  a  sinusoidal  input  signal  from 
a  Hewlett-Packard  Model  200C  audio  oscillator • 

Calibration  of  these  components  was  conducted  before 
and  after  data  takings  no  time  variation  of  response  was 
observed.  The  initial  calibration  was  performed  for  both 
instruments  by  measuring  both  input  and  output  signals  with 
the  Ballantine  VTVM.  This  necessitated  removing  the  volt- 
meter from  the  input  circuit  to  measure  output  voltages  from 
the  component  being  calibrated.  While  this  procedure  is  not 
considered  sound  practice,  the  results  of  the  initial  cali- 
bration were  corroborated  by  subsequent  calibrations „  Th® 
low  frequency  limit  of  the  audio  oscillator  used  as  a  signal 
generator  is  16  cps.  Hence,  the  equipment  could  not  be 
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calibrated  for  frequencies  below  this  limit.  The  initial 
calibrations  were  conducted  from  16  cps  in  harmonic  increments 
of  60  cps  to  an  upper  limit  of  720  cps  with  continuous  oscil- 
loscope monitoring.   Gain  factors  were  computed  for  the  pre- 

» 

amplifier,  and  for  the  amplifier  at  gain  settings  of  20,  30, 
50,  70,  and  100. 

The  final  calibrations  were  conducted  using  the  Ballan- 
tine  VTVM  for  input  measurements  and  the  thermocouple  volt- 
meter for  output  measurements.  Simultaneous  measurements  were 
taken  without  disturbing  the  calibration  circuitry.  The  two 
meters  were  previously  calibrated  against  each  other  with  a 
sinusoidal  input  signal  over  the  entire  frequency  range. 
Finally,  the  preamplifier  and  amplifier  were  calibrated  to- 
gether as  a  system  from  16  cps  to  £000  CpS  with  amplifier 
gain  settings  of  20  and  30  (i.e.,  the  gain  settings  actually 
used  during  data  taking). 

The  results  of  these  calibrations  are  presented  in  Fig. 
A-4  as  a  plot  of  relative  attenuation  or  loss  in  gain  in  db 
versus  frequency. 

The  measured  gain  of  the  preamplifier  is  seen  to  exhibit 
some  attenuation  above  35°  epsi  the  response  of  the  amplifier 
is  essentially  flat  up  to  75>0  cps,  which  is  the  upper  limit 
of  the  Type  762-B  Vibration  Analyzer.  The  preamplifier  gain 
was  taken  as  constant  over  the  operating  range  of  the  in- 
vestigation at  a  value  G  =  O.i^O,  which  was  the  arithmetic 
mean  value.  The  amplifier  gain  was  G  =  20  with  a  gain 
setting  of  20  and  Ga  =  30  with  a  gain  setting  of  30. 
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FIGURE  A-4  Component  and  System  Gain  Response 
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The  electrometer  preamplifier  was  calibrated  and 
operated  throughout  with  a  screen  grid  bias  of  +2„0  volts. 
A  circuit  diagram  of  the  electrometer  preamplifier  is  pre- 
sented in  Fig.  A-5. 
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FIGURE  A-5  Electrometer  Preamplifier  Circuit  Diagram 
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APPENDIX  B 

VELOCITY,  TURBULENCE  INTENSITY  AND 

PRESSURE  DISTRIBUTION  PROFILES 
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FIGURE  B-2  Wake  Structure  Summary  Plot,  Plate  B 
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FIGURE  B-23  Wake  Symmetry  Plot,  Plate  A 
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FIGURE  B-24  ¥ake  Symmetry  Plot,  Plate  B 
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FIGURE  B-25     Wake   Symmetry  Plot,   Plate   C 
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APPENDIX  C 
NORMALIZED  ENERGY  DENSITY  SPECTRA 
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FIGURE  0-1   Energy  Density  Spectra  at  y/t  =  0,  Plate  A 
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FIGURE  C-2  Energy  Density  Spectra  at  y/t  =  0,  Plate  B 
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FIGURE  0-3     Energy  Density  Spectra  at  y/t  =  0,   Plate   0 
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FIGURE  0-4     Energy  Density  Spectra  at  y/t  =0.5,   Plate  A 


10 


-2 


10 


-3 


E(n) 


10 


-4 


10-5 


-» ■■ 

0 


•    00 


■•,► 


t 


i* 


>4i* 


*  •  °°  ft  ° 

•    *•«.    o 


O  x/t 

•  x/t 

A  x/t 

▲  x/t 


=  9.75 
=  6.50 
=  3.55 
=   1.50 


j i    i   i  i  i  i 


6 


«  o<\> 


-**- 


1-? 


lV 


J I I    I    I  I  I 


A 

1 
o 


j I I I i    i    |   i  t 


1  11 


10 


100 


1000 


FIGURE  0-5     Energy  Density  Spectra  at  y/t  =0.5,    Plate   B 
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FIGURE  0-6  Energy  Density  Spectra  at  y/t  =  o.5,  Plate  C 
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FIGURE  0-7  Energy  Density  Spectra  at  y/t  =  0,  x/t  =  1 . 50 
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FIGURE  0-8  Energy  Density  Spectra  at  y/t  =  0,  x/t  =  3.25 
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FIGURE  C-9  Energy  Density  Spectra  at  y/t  =  0,  x/t  =6.50 
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FIGURE  0-10  Energy  Density  Spectra  at  y/t  =  0,  x/t  =  9.75 
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FIGURE  0-11   Energy  Density  Spectra,  y/t  =  0.5,  x/t  =  1.50 
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FIGURE  0-12  Energy  Density  Spectra,  y/t  =  0.5,  x/t  =  3.25 
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FIGURE  0-13  Energy  Density  Spectra,  y/t  =  0. 5,  x/t  =  6.50 
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FIGURE  0-14  Energy  Density  Spectra,  y/t  =  0.5,  x/t  =  9.75 
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FIGURE  0-15  Energy  Density  Spectrum,  Empty  Test  Section 
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APPENDIX  D 


HELMHOLTZ  RESONANCE  EFFECTS 
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1.   Plate  Spindle  Cavities. 

During  measurements  incident  to  determination  of  the 
characteristics  of  the  ambient  tunnel  turbulence  spectrum 
a  possible  Helmholtz  resonance  effect  was  observed.  With 
the  desired  flow  established,  the  amplified  turbulence  in- 
tensity voltage  was  observed  to  be  U.0  mv  rms.  Monitoring 
the  entire  available  spectrum  with  the  vibration  analyzer 
showed  no  significant  peaks.  The  test  plate  had  been  re- 
moved from  the  tunnel  test  section;  no  plugs  had  been  placed 
in  either  of  the  spindle  cavities. 

After  several  minutes  the  turbulence  voltage  level  began 
to  rise,  reaching  approximately  90  mv  rms,  and  a  discrete 
quantity  of  energy  was  observed  with  the  vibration  analyzer 
at  about  k.0   cps.  The  remainder  of  the  spectrum  indicated  an 
essentially  constant  level,  down  about  6  db  from  the  peak 
reading.   The  rms  voltage  level  gradually  increased  and  the 
peak  frequency  was  gradually  noted  to  decrease,  reaching  as 
low  as  33  cps. 

At  this  point  a  large  accumulation  of  air  bubbles  was 
noted  in  the  upper  spindle  cavity  under  the  top  cover  plate. 
The  vent  was  opened,  and  as  the  bubbles  vented  off,  the  rms 
voltage  level  decreased  to  about  1+0  mv  and  the  discrete  peak 
disappeared. 

This  phenomenon  indicated  the  possibility  of  Helmholtz 
resonance  to  the  authors.   Clearly  only  the  air  bubbles  were 
causing  the  low  frequency  energy,  but  even  a  water-filled 
cavity  might  be  expected  to  resonate  and  contribute  energy  to 
the  flow.   Even  with  the  test  plate  installed  and  the  mounting 
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spindles  passing  through  the  upper  and  lower  tunnel  walls, 
cavities  of  complicated  geometry  exist.   See  Pig.  D-l  below, 


FIGURE  D-l.  Sketch  of  Spindle  Cavity 


The  fundamental  resonant  frequency  of  a  Helmholtz  type 

r-a-ai 
cavity  can  be  estimated tJJJ  using  the  electrical  analogy  of 

a  simple  tank  circuit.  For  a  cavity  whose  dimensions  are 

known: 

S  =  surface  area  of  the  cavity  opening, 

V  =  volume  of  the  cavity, 

le  =  effective  length  of  the  "neck"  of  the  cavity, 

approximated  by  actual  length  +  0.8,y^,  and 

c  =  sonic  velocity  in  the  medium,  taken  as  1+79U 

fps  in  fresh  water,  uncorrected  for  rigidity 

of  cavity  "neck. " 
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Then  the  Hplmholtz  resonant  frequency,  nfl»  is  approximated  by: 

and  in  this  case  from  actual  measurements: 

It  is  well  known  that  the  theoretically  derived 
frequency  is  an  upper  limit  and  that  in  general  any  number 
of  effects  can  reduce  the  derived  frequency.  Furthermore, 
Morse  points  outl-,JJ  that  the  lumped  circuit  approximation 
is  less  valid  for  acoustical  systems  than  it  is  for  electrical 
circuits  because  of  the  shortness  of  the  wave  lengths  in- 
volved.  He  notes  that  Eq.  (D-l)  is  of  questionable  validity 
for  frequencies  above  1,000  cps,  for  when  a  system  is  longer 
than  a  half  wave  length  it  begins  to  behave  like  a  trans- 
mission line  and  the  wave  equations  become  involved. 
2.  The  Turbulence  Probe  Tip  as  a  Helmholtz  Resonator. 

The  above  discussion  led  naturally  to  a  consideration  of 
the  turbulence  probe  tip  as  a  Helmholtz  resonator.  From  the 
interior  dimensions  of  the  probe  tip  as  cited  in  Reference  (13) • 
and  using  Eq.  (D-l),  it  was  found  that  for  the  turbulence 
probe  tip: 

nH  ~  7060  cps 

It  is  considered  probable  that  signals  above  this  fre- 
quency will  be  attenuated  at  the  rate  of  approximately  20  db 
per  decade. 
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APPENDIX  E 
REPRESENTATIVE  SAMPLES  OP  DATA 
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SAMPLE  WAKE  SURVEY  DATA 


PLATE 

C.    x/t  =  9.75 

DATE:    3-28-61 

AMPLIFIER  GAIN:    30 

CONTRACTION  MANOMETER:    i 

*hM  =  21.30 

■    AhHg  : 

=   1.27 

VALVE 

INDICATORS   - 

FOOT 

VALVE : 

1»9 

PUMP 

DISCHARGES 

22.5 

CIRCUIT   STOP: 

12.5 

TURBULENCE  VOLTAGE    (RMS) 

TOTAL  : 

BEAD 

STATIC 

HEAD 

Thermocouple 

5     Ballantine      Pressure 

Atmos. 

Pressure 

Atmos. 

Side 

Side 

Side 

Side 

y/t 

mv 

mv 

"Hg 

"Hg 
25.86 

"Hg 

"Hg 

+  1.5 

110 

93 

20.73 

20.55 

24.47 

+1.4 

115 

102 

20.74 

25.85 

20.56 

24.47 

+1.3 

120 

105 

20.75 

25.84 

20.56 

24.46 

+1.2 

125 

110 

20.76 

25.84 

20.56 

24.46 

+1.1 

135 

122 

20.78 

25.82 

20.57 

24.45 

+1.0 

139 

127 

20.80 

25.77 

20.57 

24.45 

+0.9 

145 

130 

20.82 

25.75 

20.57 

24.45 

+0.8 

151 

132 

20.83 

25.74 

20.57 

24.45 

+0.7 

153 

136 

20.84 

25.73 

20.58 

24.44 

+0.6 

156 

138 

20.86 

25.72 

20.58 

24.44 

+0.5 

157 

138 

20.87 

25.71 

20.58 

24.44 

+0.4 

155 

142 

20.88 

25.70 

20.58 

24.44 

+0.3 

155 

141 

20.90 

25.68 

20.59 

24.43 

+0.2 

157 

143 

20.90 

25.68 

20.59 

24.43 

+0.1 

154 

142 

20.90 

25.68 

20.59 

24.43 

0 

150 

138 

20.90 

25.68 

20.60 

24.42 

-0.1 

155 

143 

20.90 

25.68 

20.59 

24.43 

-0.2 

157 

145 

20.88 

25.70 

20.59 

24.43 

-0.3 

154 

140 

20.86 

25.73 

20.59 

24.43 

-0.4 

149 

142 

20.86 

25.73 

20.59 

24.43 

-0.5 

149 

139 

20 .  82 

25.76 

20.57 

24.45 

-0.6 

150 

138 

20.80 

25.78 

20.56 

24.46 

-0.7 

147 

132 

20.80 

25.78 

20.56 

24.46 

-0.8 

142 

130 

20.77 

25.81 

20.56 

24.47 

-0.9 

135 

123 

20.76 

25.83 

20.56 

24.47 

-1.0 

133 

123 

20.75 

25.84 

20.56 

24.47 

-1.1 

122 

110 

20.74 

25.85 

20.55 

24.47 

-1.2 

117 

103 

20.74 

25.85 

20.55 

24.47 

-1.3 

115 

95 

20.73 

25.86 

20.55 

24.47 

-1.4 

100 

87 

20.72 

25.86 

20.55 

24.47 

-1.5 

97 

82 

20.72 

25.87 

20.56 

24.48 
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SAMPLE  ENERGY  SPECTRUM  DATA 


PLATE  Cf 

x/t  =  3.25. 

y/t  =0.5 

DATE:  3-29-61 

AMPLIFIES 

1   GAIN:  20 

AVERAGE  RMS  TURBULENCE  VOLTAGE 

-  THERMOCOUPL 

E:  131 

mv 

BALLANTINE 

115 

mv 

ANALYZER 

SENSITIVITY: 

:  10 

ANALYZER 

SELECTIVITY: 

:  SHARP 

ANALYZER 

CALIBRATION 

CONSTANT: 

5.31  x  105 

vo 

Its"1 

CONTRACT] 

!0N  MANOMETER:  Ahu  =  21 

.35o  Ahrr>.  = 

=  1 

,36 

•50 
-  20. 

A*itotal  = 

:  25.62  -  20, 

■95     ^static  =  2^ 

66 

FREQUENCY  ANALYZER 

FREQUENCY 

ANALYZER 

FREQUENCY  ANALYZER 

(cps) 

READING 

(cps) 

READING 

(cps) 

READING 

n 

Rd 

n 

Rd 

n 

Rd 

2.5 

7.5 

70 

22 

230 

28 

3 

9.4 

75 

25 

240 

26 

4 

8.9 

80 

28 

250 

24 

5 

9.5 

35 

31 

260 

25 

6 

11.0 

90 

38 

270 

24 

7 

9.9 

95 

60 

280 

27 

8 

11.7 

98 

75 

290 

26 

9 

12.1 

100 

92 

300 

27 

10 

12.5 

102 

95 

320 

25 

11 

12.6 

104 

88 

340 

25 

12 

12.6 

106 

75 

360 

25 

13 

12.4 

108 

70 

380 

25 

14 

12.4 

110 

60 

400 

25 

15 

12.5 

115 

45 

420 

24 

17 

13.4 

120 

39 

440 

23 

20 

11.8 

130 

30 

460 

22 

22 

12.2 

140 

28 

480 

22 

25 

12.8 

150 

26 

500 

22 

28 

14.0 

160 

27 

520 

22 

30 

14.8 

170 

27 

540 

22 

35 

14.8 

180 

26 

560 

22 

40 

17.0 

190 

28 

580 

21 

45 

17.8 

195 

30 

600 

21 

50 

19.6 

200 

33 

650 

20 

55 

20 

205 

32 

700 

20 

60 

21 

210 

31 

750 

20 

65 

20 

220 

30 

NOTE:  Analyzer  readings  from  n  =  2.5  through  n  =  50  are  time 

averages  of  a  number  of  individual  readings  taken  at  fixed 
time  increments  in  an  attempt  to  smooth  out  data  scatter. 
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